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19.  (cont.) 

The  NORSAR  Detection  Processing  system  has  been  operated  through  the 
reporting  period  with  an  average  uptime  of  99.6  per  cent.  A  total  of 
2028  seismic  events  have  been  reported  by  NORSAR  in  the  period.  The 
performance  of  the  continuous  alarm  system  and  the  automatic  bulletin 
transfer  by  telex  to  AFTAC  have  been  satisfactory.  Processing  of 
requests  for  full  NORSAR/NORESS  data  on  magnetic  tape  has  progressed 
without  problems. 

Investigations  into  further  potential  improvements  in  the  NORSAR 
array  processing  system  have  continued.  A  new  Detection  Processor  (DP) 
program  has  been  developed  and  tested  in  an  off-line  mode.  This 
program  is  flexible  enough  to  conduct  both  NORSAR  and  NORESS  detection 
processing  as  is  done  today,  besides  incorporating  improved 
algorithms.  The  program  will  be  implemented  after  some  further  testing 
and  refinements  have  been  made. 

Field  maintenance  activity  has  Included  regular  preventive  maintenance 
at  all  subarrays  and  occasional  corrective  actions  as  required.  No 
special  problems  have  been  noted  in  the  performance  of  the  field 
installations.  In  addition,  field  maintenance  personnel  have  par¬ 
ticipated  in  experiments  involving  the  NORESS  array  as  well  as  a 
planned  second  small  aperture  array. 

A  wide-band  slowness  estimation  technique  has  been  investigated  by 
processing  data  from  several  events  from  the  same  location.  Ten  quarry 
blasts  at  a  dam  construction  site  in  western  norway ,  ten  explosions 
from  a  mining  area  in  western  Russia  and  sixteen  Semipalatinsk  nuclear 
explosions  were  selected.  The  major  conclusion  from  this  study  is  that 
employing  a  wider  frequency  band  clearly  tends  to  increase  the  stabi¬ 
lity  of  the  slowness  estimates,  provided  the  signal-to~noise  ratio  is 
adequate  over  the  band  of  Interest.  The  stability  yas  found,  par¬ 
ticularly  for  Pn,  to  be  remarkably  good  for  the  western  Norway  quarry 
blasts  when  using  a  fixed  frequency  band  for  each  phase  for  all  ten 
events. iXThis  approach  also  contributed  a  full  separation  between  Sn 
and  Lg  phases  in  slowness  space,  with  associated  average  phase  veloci¬ 
ties  of  4.71  km/s  and  3.91  km/s,  respectively.  This  clear  separation 
could  not  be  obtained  for  the  western  USSR  events,  part'y  due  to 
narrower  usable  bandwidth.  On  the  other  hand,  the  Semipalatinsk  explo¬ 
sions,  at  teleseismic  distances,  showed  very  consistent  results  when 
applying  the  wide-band  technique. 
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The  regular,  automatic  recording  of  noise  spectra  has  now  been  taking 
place  for  close  to  one  year.  Both  NORESS  and  High  Frequency  Element 
data  have  been  processed  with  a  view  to  obtaining  quantitative  esti¬ 
mates  of  noise  level  variations  on  a  seasonal  and  diurnal  basis.  Main 
results  are: 

a)  The  average  noise  spectrum  falls  off  rapidly  (~  f"^  or  24  d B/ 
octave)  between  0.5  and  2  Hz.  Above  2  Hz,  and  up  to  at  least 
50  Hz,  the  spectral  falloff  is  ~  f-^  (or  12  dB/octave). 

b)  Over  a  one-year  period,  the  standard  deviation  of  the  noise  level 
(measured  in  dB)  is  relatively  constant  across  the  entire  fre¬ 
quency  band  considered  (0.5-50  Hz),  and  is  typically  3-4  dB.  There 
is,  however,  a  tendency  toward  increased  noise  standard  deviation 
both  at  the  low  end  and  at  the  high  end  of  this  band. 

c)  On  the  average,  there  is  a  systematic  diurnal  variation  in  noise 
level,  particularly  on  workdays.  The  typical  difference  between 
nighttime  and  daytime  noise  (measured  at  02  and  13  hours  local 
time,  respectively)  is  5  dB  in  the  band  2-50  Hz.  At  frequencies 
below  1.5-2  Hz,  the  difference  is  insignificant.  At  particular 
frequencies,  e.g. ,  6  Hz,  there  are  peaks  in  the  noise  spectra 
caused  by  specific  industrial  activity,  ana  the  corresponding 
diurnal  variation  may  at  the  extreme  exceed  10  dB. 


A  comprehensive  analysis  has  been  carried  out  using  data  recorded  by 
an  experimental  small-aperture  array  (F1NESA)  in  Finland,  and  com¬ 
paring  the  results  to  those  obtained  from  NORESS.  We  have  found  that 
the  noise  correlation  properties  closely  resemble  those  found  at  the 
NORESS  site;  thus  the  design  philosophy  used  for  NORESS  in  order  to 
exploit  negative  noise  correlation  seems  to  have  more  general 
applicability.  Otherwise,  the  noise  level  at  the  Finnish  site  (north 
of  Helsinki)  is  comparable  to  that  of  NORESS  at  1  Hz,  and  slightly 
higher  than  for  NORESS  at  higher  frequencies.  .Joint  processing  of 
regional  events  recorded  by  both  FINESA  and  NORESS  has  shown  that 
there  is  a  potential  for  significant  improvement  in  event  location 
using  two  arrays  rather  than  one.  The  RONAPP  detection  algorithm  has 
been  applied  to  FT.NESA  data  for  a  two-week  period,  and  has  produced 
quite  good  results.  In  particular,  it  has  been  possible  to  obtain 
automatic  separation  of  P  phases  and  secondary  phases  with  close  to 
100  per  cent  accuracy. 

A  recursive  method  to  model  multiple  scattering  by  topography  on  an 
interface  or  surface  has  been  developed.  The  need  for  such  a  method 
has  arisen  especially  in  situations  where  the  scale  of  lateral 
variations  is  of  the  order  of  a  wavelength.  At  this  stage  the  method 
has  been  coded  for  a  solid-liquid  interface.  Initial  results  have 
demonstrated  a.o.  the  range  of  applicability  ot  first-order 
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approximations.  Having  in  mind  the  observed  frequency  dependent  wave- 
number  anomalies  at  NORESS,  we  plan  to  apply  the  method  to  scattering 
by  surface  topography  as  well.  This  will  enable  us  in  the  future  to 
compare  the  frequency  dependent  phase  perturbations  to  the  usual  sta¬ 
tic  corrections  due  to  topography. 

A  large-scale  study  has  been  completed  with  the  aim  to  determine  phase 
characteristics  and  detectability  of  mine  explosions  and  rock  bursts 
in  Polish  mining  areas.  We  have  obtained  a  data  base  of  more  than  500 
extremely  well-located  seismic  events  (hypocentral  uncertainty  a  few 
tens  of  meters)  of  magnitude  range  1.5-3. 5  (Ml).  Conclusions  from  the 
study  are  that  RONAPP  automatic  P-detection  is  excellent  (90  per  cent 
threshold  Ml  2. 7-2. 8  for  the  Lubin  area  at  1060  km  distance),  whereas 
automatic  detection  of  secondary  phases  seldom  occurs.  The  main  reason 
is  that  the  Tornquist-Teisseyre  tectonic  zone  in  Poland  acts  as  a 
barrier  to  Lg  propagation.  Using  NORESS-recorded  P  and  S  phases  in 
conjunction  with  S  azimuths  gives  reasonable  location  estimates  for 
those  events  where  both  phases  can  be  seen;  typically  the  location 
error  is  50-60  km. 

The  accuracy  of  estimated  arrival  times  and  azimuths  depends  strongly 
on  the  shape  and  corner  frequency  of  recorded  phases.  Characteristic 
patterns  of  P-wave  records  might  provide  important  information  into  an 
expert  system  approach  to  regional  array  processing. 
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I .  SUMMARY 

This  Semiannual  Technical  Summary  describes  the  operation,  maintenance 
and  research  activities  at  the  Norwegian  Seismic  Array  (NOKSAR)  for 
the  period  1  Oct  1986  -  31  Mar  1987. 

The  uptime  of  the  NORSAR  online  detection  processor  system  has 
averaged  99.6  per  cent  during  the  reporting  period  as  compared  to  99.3 
for  the  previous  period.  There  have  been  a  number  of  problems  with  the 
communication  lines  to  subarray  04C  for  most  of  the  period.  Otherwise 
the  system  has  generally  shown  a  reliable  performance.  A  total  of 
2028  seismic  events  were  reported  in  the  NORSAR  monthly  seismic  bulle¬ 
tins  in  this  period,  giving  a  daily  average  of  11.1  events. 

The  new  Detection  Processor  (DP)  program  has  now  been  tested  and  found 
to  perform  satisfactorily.  This  DP  system  is  described  in  some  detail 
and  some  preliminary  results  from  processing  done  so  far  are 
documented . 


Field  maintenance  activity  has  Included  regular  preventive  maintenance 
at  all  subarrays  and  occasional  corrective  actions  as  required.  No 
special  problems  have  been  noted  in  the  performance  of  the  field 
installations.  In  addition,  field  maintenance  personnel  have  par¬ 
ticipated  In  planning  for  the  installation  of  the  NORESS  11  array  in 
Finnmark. 


The  research  activity  is  summarized  in  section  VII.  Section  VII, l 
describes  wide-band  slowness  estimation  using  a  small  aperture  seismic 
array.  Section  VII. 2  presents  NORESS  noise  spectral  studies,  with  spe¬ 
cial  emphasis  on  noise  level  characteristics.  In  Section  VII. 3  initial 
results  from  data  analysis  using  the  FINESA  experimental  small  aper¬ 
ture  array  are  presented.  Section  Vtt.4  describes  a  method  to  account 
for  multiple  scattering  due  to  topography  of  internal  boundaries  or  of 
the  free  surface.  In  Section  VII. 5  NORESS  capability  for  detection  and 
location  of  mining  tremors  in  the  Lubin  area  in  Poland  is  discussed. 
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II.  OPERATION  OF  ALL  SYSTEMS 

II.  I _ Detection  Processor  (DP)  Operation 

There  have  been  90  breaks  in  the  otherwise  continuous  operation  of  the 
NORSAK  online  system  within  the  current  6-month  reporting  interval. 

The  uptime  percentage  for  the  period  is  99.6  as  compared  to  99.3  for 
the  previous  period. 

Fig.  It. 1.1  and  the  accompanying  Table  It. 1.1  both  show  the  daily  DP 
downtime  for  the  days  between  1  October  198b  and  31  March  1987.  The 
monthly  recording  times  and  percentages  are  given  in  Table  11.1.2. 


The 

breaks  can  be  grouped  as 

follows 

a) 

Hardware  failure 

b) 

Stops  related  to  program 

work  or 

c) 

Hardware  maintenance  stops 

d) 

Power  jumps  and  breaks 

e) 

TOD  error  correction 

O 

Communication  lines 

The  total  downtime  for  the  period  was 
mean-time-between-failures  (MTBF)  was 
the  previous  period. 


1 

error  0 
3 
3 

12 

39 

16  huurs  and  14  minutes.  The 
2.0  days,  as  compared  to  2.3  for 


J.  Torstveit 


Month 


Month 

DP  uptime 

hours 

DP  uptime 

X 

No .  of 

DP  breaks 

No.  of  days 

with  breaks 

DP  MTBF* 

(days) 

OCT 

737.33 

99.1 

13 

12 

2.2 

NOV 

713.78 

99.4 

19 

14 

1.6 

DFC 

743.73 

99.9 

13 

9 

2.2 

JAN 

742.98 

99.9 

21 

7 

1.4 

FEB 

668.26 

99.4 

7 

4 

3.5 

MAR 

743.65 

99.9 

18 

9 

1.6 

99.6 

90 _ 

35 

2.1 

*  Mean-time-between-failures  »  total  uptime/no.  of  up  Intervals. 


Table  U.1.2  Online  system  performance,  1  Oct  1986  -  31  Mar  1987 
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Table  II.  1. 1  Dally  DF  downtime  in  the  period  1  Oct  I9bb  -  31  Mar  19«7. 


Table  tt.1.1  (cont 
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Array  communication;: 


Reference  is  made  to  Table  VII.2.1  in  connection  with  communication 
system  performance  in  the  period. 

The  most  affected  systems  have  been  04C  and  06C,  04C  due  to  reduced 
line  performance  and  06C  due  to  an  assumed  intermittent  operating 
equalizer. 

Besides,  a  nonlocated  noise  source  has  affected  all  communications 
systems,  especially  during  night  hours.  Although  not  causing  an 
excessive  error  rate,  the  irregularity  caused  many  resynchronizations 
of  the  communications  lines.  Xt  is  clear  that  the  noise  is  created 
somewhere  along  the  NTA's  carrier  systems.  All  our  communications 
systems  are  located  within  the  same  channel  group.  On  27  April,  NT A 
switched  to  a  different  equipment  group,  and  this  particular  problem 
then  disappeared. 

Otherwise,  and  as  the  table  indicates,  most  systems  have  been 
reliable. 

NTA/Llllestrom  and  Haraar  wilJ  take  a  new  initiative  medio  April 
regarding  04C.  The  line  will  be  remeasured  and  if  a  certain  degree  of 
improvement  in  the  line  condition  is  found,  new  modems  will  be  hooked 
up  to  see  if  their  equalizers  can  handle  the  line  distortion.  The 
modems  will  then  be  kept  until  the  line  is  repaired. 

Summary 

Oct  86  03C  had  been  unreliable  for  some  time,  continued  week  40. 

Apart  from  an  irregularity  week  42,  the  03C  communication 
system  remained  stable  throughout  the  period. 
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02B  was  affected  week  43  (19. 7%)  caused  by  an  assumed 
line  fault  toward  the  subarray. 

The  performance  of  the  remaining  systems  was  acceptable. 

Average  communication  system  outage/degradation  5  weeks 
all  systems:  1.67%  and  0.11%  (-  020). 

Nov  86  Apart  from  04C,  reliable  performance  all  systems.  Although 
having  a  significant  amount  of  errors  (0.93%)  subarray  04C 
was  not  permanently  masked. 

Average  communication  system  outagi/degradation  4  weeks 
all  systems:  0.13%  and  0.02%  (-  040). 

Dec  86  04C  still  out  of  operation  due  to  changed  cable/line 

characteristics.  NTA/Hamar  visited  the  subarray  2  Dec.  Level 
measurements  revealed  too  high  input  level  (-13  dBm),  which, 
unless  otner  Irregularities  also  existed,  isolated  should  not 
cause  problems. 

NMC  visited  the  subarray  12  and  18  Dec.  Attempts  were  made  to 
replace  n oderas  and  loop  control  lo;ic,  but  without 
improvements.  Communications  lines  between  Kjeller  and  the 
subarray  were  measured  with  respect  to  attenuation  distortion 
for  the  actual  frequencies  300-2400  Hz.  Between  300  and 
800  Hz  the  attenuation  was  0.6  and  1.0  dB  m  in  the  directions 
Kjeller-ITV  and  CTV-Kjeller,  respectively.  Between  800  and 
2400  Hz  the  attenuation  CTV-Kjeller  was  0.5  dBm,  while  be¬ 
tween  Kjeller  and  the  CTV  the  attenuation  was  2.2  dBm,  the 
latter  normal. 

Late  December  NTA/Elverum  swapped  cable  pairs  between 
Berg  Central  and  the  CTV  without  improving  the  quality. 

Due  to  very  low  temperatures  and  < eep  snow  in  the  subarray 
area,  no  further  actions  were  expi cted  until  conditions 
changed. 
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060  was  affected  frequently  week  52,  causing  high  outage 
figures.  An  inr ;rmi r tent  operating  equalizer  was  probably  the 
reason.  After  29  De  ember  the  performance  ha;  been  reliable. 

Average  communication  system  outage/degradat Lon  A  weeks 
ail  systems:  11. 62,  and  0.02%  (-04C,  06C). 

Jan  87  It  has  been  necessary  to  mask  01B  and  03C  a  few  times  due  to 
spikes  in  data  (between  5  and  12  Jan). 

04C.  No  progress  wi'h  regard  to  the  communication  system 
outage.  NTA/Hamar  has  checked  equalizers  (Elverum),  and 
cables  between  Hernts  and  Berg  Central.  Cables  between  Berg 
Central  and  the  subarray  remain  to  be  measured  thoroughly. 

02C  was  affected  27  January  in  connection  with  planned  cable 
work  in  Mesnalien  (NTA/Hamar). 

In  connection  with  expansion/ improvement  of  the  telecom¬ 
munication  circuits  in  the  Kjeller-Lillestr0m  area  most 
institutions  were  affected,  including  NORSAR  and  the  neigh¬ 
boring  institution  RBK. 

In  our  operator/modera  office,  the  data  and  telecommunication 
junction  box  was  replaced  and  cables  were  cut /spliced.  A  new 
cable  (20  pairs)  was  installed  between  the  junction  box  and 
RBK  telephone  central.  The  above-mentioned  activities  took 
place  between  28  January  and  1  February,  and  all  our  data  and 
telecommunication  circuits  were  broken  for  short  periods. 

06C.  Communication  circuit  instability  contirued  through  2U 
December,  but  reliable  performance  after  2  January. 

Average  communication  system  outage/degradation  5  weeks  all 
systems  (-  04C):  0.17%. 
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Feb  87  04C  still  out  of  operation  as  NTA/l.amar  was  not  able  to 

repair  th«  cable  between  Berg  Central  and  the  CTV. 

NMC  visit* d  the  subarray  27  February  and  several  tests  were 
carried  out  in  cooperation  with  NDI'C.  Error  rate,  especially 
in  the  di  ection  towards  Kjeller,  was  high. 

Average  communication  system  outag-  /degradation  4  weeks: 

0.02%  (-  i4C). 

Mar  87  NTA/Hamar  and  Lillestr^m  ready  for  a  new  initiative  mid- 

April  in  onnection  with  04C.  The  omplete  communication  path 
will  be  r  -measured,  tf  positive  ch  inges  in  line  quality  since 
last  meas  irements ,  new  modems  (wit  i  automatic  equalization) 
will  be  h  >oked  up  in  the  CTV  and  h  :re  at  the  NDPC.  If  usable 
results  a  e  obtained  with  the  mode  as ,  they  will  be  kept  until 
cables  or  other  parts  of  the  commu  ilcation  system  causing  the 
trouble  a  e  repaired. 

06C.  Err >r  figures  for  March  (13.80%)  were  also  very  high, 
mainly  du  :  to  frequent  resynchroni '.at ion  of  the  communication 
system. 

Since  23  larch  all  NORSAK  communic rtlon  systems  were  affected 
by  a  nonl  icated  noise  source,  mainly  during  night  hours. 
Although  aoderate  error  rate,  such  conditions  cause  unwanted 
line  res) achronization,  and  extra  load  on  the  communication 
processor.  This  problem  persisted  until  27  April,  when  NT  A 
switched  their  channel  group  equipnent  for  all  lines  to  a 
spare  channel  group. 

Average  c  rmraunication  system  outape/degradation  3  weeks 
(10-12),  all  systems  (-  04C):  0.13  %.  4  v:eeks  (10-13), 

01A-03C  •-  04C ,  06C) :  0.77%. 


O.A.  Hansen 
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Ttt.  ARRAY  PERFORMANCE 


Event  Processor  Operation 

In  Table  III# I  some  monthly  statistics  of  tne  Event  Processor  opera¬ 
tion  are  given: 


Teleseismic 

Core  Phases 

Sum 

Daily 

OCT  86 

195 

255 

450 

14.5 

NOV  86 

224 

126 

350 

11.7 

DEC  86 

175 

145 

320 

10.3 

JAN  87 

205 

74 

279 

9.0 

FEB  87 

188 

107 

295 

10.5 

MAR  87 

244 

90 

334 

10.8 

1231 

797 

2028 

11.1 

Table  tll.l  Event  Processor  Statistics,  Oc'.ober  198b  -  March  1987. 


8.  Piulsen 
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IV.  IMPROVEMENTS  AND  MODIFICATIONS 

In  NORSAR  Scientific  Report  Mo.  1-86/87  we  outlined  an  implementation 
plan  for  the  upgrade  of  the  NORSAR  processing  system. 

Part  1,  to  implement  new  NORSAR  offline  detection  system,  has  been 
completed,  and  tested  parallel  to  the  present  recording  and  detection 
system.  This  new  program  system,  DP  -  Detection  Processor  -  has  been 
found  to  perform  most  satisfactorily  the  NORSAR  detection  task. 

In  the  'old'  system,  the  detection  processing  is  divided  into  two 
steps.  A  MODCOMP  computer  system  communicates  with  the  NORSAR 
subarrays  and  performs  filtering  and  subarray  beamiorraing  on  the  data. 
Thereafter,  the  data  is  transmitted  each  0.5  second  to  the  IBM 
4341/L01  processor,  where  ariay  beaoforming  and  STA/LTA  detection  is 
done.  The  new  system  perforus  all  operations  in  one  program,  to  be 
run  on  the  IBM  4341.  As  opposed  to  the  'old'  system  where  array  beam- 
forming  is  done  using  subarray  beams,  the  new  system  forms  beams  using 
individual  delays  for  all  tht  42  instruments.  Incoherent  beamforming 
is,  however  done  using  subarray  beams. 

In  this  report  we  will  describe  some  of  the  details  of  the  DP  system 
and  document  some  preliminary  results  from  the  processing  done  so  far. 

In  Scientific  Rep.  No  1-86/87,  we  outlined  some  general  rules  for 
developing  processing  programs  and  gave  examples  of  application 
of  two  parts  of  a  data  analysis  program:  the  data  retrieval  part  and 
the  presentation  of  the  data  or  the  results.  The  actual  analysis  part 
in  between  these  steps  consists  of  applying  a  variety  of  processing 
techniques,  and  we  will  classify  these  into  two  main  categories: 
"segmented"  and  "continuous". 
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We  define  segmented  processing  as  applying  some  analysis  to  a  segment 
of  data  and  present  results  for  that  segment.  Continuous  processing  is 
defined  as  a  repeatable,  cyclic  process  applied  to  a  data  stream  that 
in  principle  is  of  infinite  length. 

In  a  detection  system,  the  processes  applied  (filtering,  STA/LTA 
operators,  beamforming)  need  at  any  time  point  some  of  the  'old'  data. 
Therefore,  we  cannot  produce  detector  results  by  only  processing  only 
a  few  seconds  of  the  data,  we  need  at  least  some  minutes  before  detec¬ 
tor  statistics  stabilize. 

Continuous  processing  may  be  implemented  as  segmented  processing  of 
long  data  intervals  which  follow  in  time,  aid  continuity  may  be 
secured  by  reprocessing  some  part  of  the  data.  Memory  requirements 
may  be  too  large  for  such  an  approach,  but  this  can  be  solved  by 
dividing  the  job  into  smaller  steps,  and  let  each  step  output 
Intermediate  results. 

Another  approach  is  to  process  a  smaller  se  .ment  of  data  and  retain 
some  of  the  previously  processed  data  for  continuity  when  the  next 
segment  is  brought  into  processor  memory.  We  define  this  as  "con¬ 
tinuous  segmented”  processing,  and  this  is  the  technique  implemented 
in  the  new  DP  system.  For  each  segment,  the  routines  operate  as  if 
only  one  single  segment  of  data  were  to  be  processed.  The  idea  is  that 
some  processes  then  require  only  one  prograi.  routine,  such  that 
program  parts  of  different  systems  may  use  the  same  routines. 

In  practice,  the  continuous  segmented  processing  is  implemented  by 
keeping  all  data  in  a  memory  buffer.  At  any  time  instant,  a  process 
which  uses  some  'old'  data  (e.g.  filtering  and  beamfcrming)  may 
reference  time  points  backwards.  This  is  made  possible  by  always 
keeping  a  defined  amount  of  data  as  'old'  data.  E.g.,  for  NOKSAR  pro¬ 
cessing,  the  amount  of  'old'  data  needed  at  sny  time  point  is  derived 
from  the  maximum  be.imforming  delays  and  from  the  delay  used  when  up- 
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dating  LTA.  Currently  this  parameter  is  9  seconds.  That  is,  if  6U.0 
seconds  of  data  segment  length  is  used,  then  we  always  keep  69. U 
seconds  of  data.  Each  time  a  new  segment  is  brought  in,  9.0  seconds  of 
data  is  moved.  Processing  tine  for  this  moving  of  data  is  only  about 
0.2  X  of  the  total  time  spent  on  processing  the  NORSAR  array  data. 

All  the  routines  then  process  60.0  seconds  segment  of  data. 

The  system  may  process  a  defined  number  of  segments  (steps),  or  do 
processing  as  long  as  there  ire  data.  The  DP  system  is  thus  designed 
to  allow  for  both  segmented  tnd  continuous  segmented  techniques. 

The  program  system  is  divided  into  several  FORTRAN  subroutines. 

We  have  earlier  described  the  benefits  of  defining  a  common 
internal  data  buffer  structure  for  use  at  NOkSAR.  The  new  analysis 
program  developed  by  T.  Kvaerua  (SSA  -  Seismic  Signal  Analyser),  uses 
the  same  data  read  routines  is  the  detection  system,  and  the  same 
internal  data  structure.  So  '  ar  we  may  use  all  the  data  recorded  by 
NORSAR  systems.  (NORSAR,  NORI.SS,  NORESS  High  frequency  system, 

FINESA,  Level  2  data).  CSS  data  format  will  also  be  implemented. 

Using  this  structure,  the  processing  code  is  not  tied  to  specific 
data,  and  coding  done  for  segmented  processing  in  the  SSA  system  may 
be  adopted  by  the  detection  i-ystem. 

The  principal  idea  behind  the  internal  structure  is  to  use  one 
large  buffer  for  all  data.  Instrument  data  are  called  'channels', 
likewise  any  other  vector.  The  buffer  then  contains  a  number  of  chan¬ 
nels,  each  channel  described  by  its  position  in  buffer,  number  of 
points,  name,  distance  between  points  (sample  rate)  etc. 

The  DP  system  currently  consists  of  the  following  main  routines: 
an  initialization  routine,  a  orocessor  scheduling  routine,  routine  for 
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filtering,  routine  for  beamforming,  routine  for  STA/LTA  detection, 
routine  for  incoherent  STA/LTA  detection,  and  routines  for  forming 
STA,  LTA,  and  STA/LTA  time  series.  There  are  of  course  also  routines 
for  reducing  detector  output,  and  recording  of  detections. 

The  detector  initialization  routine  will  either  read  commands  from  the 
terminal  keyboard,  or  from  a  detector  recipe  coded  by  the  user.  The 
initialization  part  will  build  up  the  data  Duffer  (by  initializing 
channel  pointers),  and  build  up  processing  orders  for  the  scheduler. 

Figs.  IV. 1,  IV. 2  and  TV. 3  show  extracts  frou  the  NORSAR  and  NORESS 
detection  recipes,  which  may  be  used  as  examples  for  the  commands  used 
by  the  current  DP  system.  The  commands  are: 

TERMINAL  -  Switch  from  disk  to  terminal  input. 

TYPE  -  Define  amount  of  screen  output  from  routines. 

STEP  -  The  number  of  segments  to  process.  Start  time 

is  given  by  main  controlling  program. 

DPTTME  -  If  start  time  for  processing  is  specified  with 

QPTTME  command,  then  the  system  will  process  the 
number  of  segments  specified,  and  at  end  of 
processing,  the  data  segment  in  the  buffers  will 
contain  that  time. 

DATA  BUFFER  LENGTH  60.0  SECONDS 

-  Set  segment  procesing  length  to  60.0  seconds. 

OLD  BUFFER  LENGTH  9.0  SECONDS 

-  Set  necessary  backup  time  to  9.U  seconds. 

SET  DETECTION  STATE  TRUE 

-  Initially  for  the  first  segment,  the  program  will  believe 
that  all  beams  have  STA/LTA  greater  than  threshold. 

SET  INITIAL  LTA  8.0 

-  Set  initial  LTA  value  of  8.0. 

SET  REFERENCE  02BO0 

-  Set  reference  station  to  02BOO. 
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SET  THRESHOLD  3.2 

-  Set  STA/LTA  threshold  to  3.2. 

SET  STX  LENGTH  0.5  SECONDS. 

-  Set  length  of  STX  to  0.5  second.  STX  is  then  the  average  of 
0.5  second  of  rectified  data.  STX's  will  be  sampled 

with  a  rate  equal  to  its  length  and  are  not  overlapping. 

SET  STA  LENGTH  1.5  SECONDS. 

-  Set  STA  length  of  1.5  second.  STA  will  be  the  average  of 
STX's  over  1.5  second.  STA  must  cover  an  Integer  number 
of  STX's  and  is  sampled  with  the  same  rate  as  STX. 

SET  LTA  SIGMA  7.0 

-  LTA  at  one  time  instance  will  be  STX* 2.0** (-SIGMA)  plus 
the  former  LTA*( l .0-2 .0**(-SlGKA) ) .  The  STX  used  is 
not  from  'current'  time  point,  but  delayed  9.0  seconds. 

(See  SET  LTA  DELAY). 

SET  LTA  BETA  b.O 

-  If  STA/LTA  is  above  threshold,  the  LTA  will  be  updated 
with  BETA  during  this  detection  state  rather  than  SIGMA. 

SET  LTA  DELAY  9.0  SECONDS 

-  See  LTA  SIGMA  comment. 

SET  REDUCTION  3.0  SECONDS 

-  Time  constant  used  for  detection  reduction.  For  an  event, 
a  great  number  of  beams  may  detect.  If  detection  start 
time  is  closer  than  3.0  seconds  to  another  detection,  then 
the  detections  are  reported  as  one  detection,  using  the 
detector  results  for  the  one  with  largest  STA/LTA. 

INPUT  THE  FOLLOWING  STATIONS: 

-  Give  a  list  of  station  names  to  get  the  data. 

FILTER  FILNAM  FILTER  FLTTYP  BANDLW  BANDHG  I0RD 

1  FBI  BU  BP  1.2  3.2  3 

-  Give  a  list  of  filters  that  may  be  used. 

Filters  are  later  referred  to  by  name  (FBI). 
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FILTER  ALL  WITH  FBI  OUTPUT  xxxxxFSl  or 

FILTER  xxx  WITH  FBI  OUTPUT  yyy 

-  Filter  all  stations  with  filter  -B1  and  give  name  to  the 
output  filtered  stations  by  usin;  the  first  5  characters 
of  the  input  (unfiltered)  station  concatenated  with  the 
filter  name.  Or  filter  station  with  name  xxx  with  filter 
FBI  (example)  and  give  output  channel  name  yyy. 

DEFINE  NAO-AB  AS  ALL  xxxxxFBl 

-  Define  a  configuration.  Give  the  name  NAO-Ab  to  that 
configuration.  In  this  example  u >e  all  channels  with 
FBI  in  the  name. 

USE  NAO-AB 

-  In  the  following,  use  configuration  NAO-AB  lor  beamforming. 

USE  BEAM-TABLE  NAO-451  TABLE  * 

-  In  the  file  :  NAO-451  TABLE  *  ,  the  program  will  find 
a  list  of  delays  for  a  number  of  beams. 

USE  MAX  AMPLITUDE  3.0 

-  If  an  input  unfiltered  instrument  has  a  maximum  absolute 
amplitude  that  is  at  least  3.0  times  the  average  of  the 
maximum  amplitudes  at  all  channels  for  a  segment,  then  treat 
that  instrument  as  inoperative  for  the  whole  segment. 

USE  NUMBER  OF  ITERATIONS  2  FOR  MAXIMUM_AMPI ITUDE  MASKING 

-  Iterate  the  maximum  amplitude  spike  detector  2  times. 
Exclude  bad  instruments  from  average  from  each  iteration. 

STATIONS 

-  Give  list  of  stations  tor  which  delays  are  given  in  a 
beam  delay  table. 

DELAYS  FOR  BEAM-TABLE  NAO-451 

-  List  of  bears  follow.  A  beam-tat  le  entry  ha;,  name  ot 
the  beam,  number  of  stations  for  which  there  are  delays, 
and  the  delays  for  those  stations.  For  one  heam-table, 
there  must  be  the  same  number  of  stations  giving  delays 
(as  defined  in  the  STATIONS  commind).  When  using  a 
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DETECT  ON  BEAM  command,  the  system  will  find  the  stations 
in  the  configuration  used,  and  the  delays  for  that 
configuration  (a  subset  of  STATIONS). 

DETECT  ON  ALL 

-  For  each  of  the  bi  ams ,  which  are  defined  so  far  by  a  USE 
BEAM-TABLE  command  ,  form  the  individual  beam,  and 

run  STA/LTA  detector  on  that  beam. 

DETECT  ON  BEAM  xxx 

-  Form  beam  xxx  and  run  STA/LTA  detector  on  that  beam, 
which  is  defined  so  far  by  a  USE  BEAM-TABLE  command. 

DETECT  ON  xxx 

-  Run  STA/LTA  detector  on  the  channel  named  xxx. 

For  incoherent  detection,  use  DETECT  INCOHERENT  ON. 

MAKE  BEAM  xxx  OUTPUT  yyy  (USING  ccc)  (VEL  v)  (AZI  a) 

-  Form  a  beam.  If  xx<  is  a  name  of  a  beam  in  a  beam-table, 
form  that  beam  using  those  delays.  Give  the  beam  name  yyy. 
Use  configuration  ccc  if  specified.  Use  given  velocity  v 
and  azimuth  a  if  specified.  If  beam  name  xxx  does  not 
exist,  give  name  x  ,x  to  a  new  beam  entry. 

MAKE  STX  OF  xxx  OUTPUT  yyy 
MAKE  STA  OF  xxx  OUTPUT  yyy 
MAKE  LTA  OF  xxx  OUTPUT  yyy 
MAKE  SNR  OF  xxx  OUTPUT  yyy 

-  Make  STX,  STA,  LTA,  or  STA/LTA  time  series  from  the 
data  channel  named  xxx.  Output  resulting  data  channel 
will  be  given  name  yyy. 

The  commands  are  used  to  build  up  processing  orders  (tasks  to  do) 
for  the  processing  scheduler,  and  must  therefore  be  given  in  the 
correct  order  of  processing  sequence.  The  user  must  decide  whether  to 
beamform  first  and  then  filter  the  beam  or  filter  before  beamformimg, 
and  do  the  things  in  correct  order.  When  initializing  the  system,  the 
program  builds  up  memory  buffers  necessary  to  handle  the  recipe,  and 
gives  error  messages  if  buffers  are  exceeded.  The  program  also  checks 
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for  consistency  with  regard  to  sample  rate,  STX/STA  integration 
lengths,  and  data  buffer  lengths. 

These  commands  are  sufficient  to  perform  all  operations  concerning  the 
STA/LTA  detector  algorithms  used  at  NORSAR.  incoherent  beamforming  may 
be  done  by  either  using  STA  or  STX  time  seties  for  beamforming  and/or 
detection.  STA  rate  may  be  made  equal  to  instrument  sample  rate  by 
setting  STX  length  to  original  sample  interval  time.  Beam  delays  are 
given  in  seconds  and  must  be  given  as  relative  arrival  times.  A  posi¬ 
tive  delay  thus  means  arrival  time  later  thin  the  reference  instru¬ 
ment.  The  program  system  itself  corrects  all  delays,  such  that  the 
instrument  with  latest  arrival  gets  zero  delay,  and  all  other  delays 
are  negative.  The  program  system  may  easily  be  adjusted  to  do  weighted 
beamforming,  do  power  STA/LTA  detection,  or  other  processes  that 
fit  a  continuous  processing  scheme. 

By  using  these  commands  we  have  programmed  i  recipe  for  doing  NORSAK 
subarray  and/or  array  beamforming,  coherent  and  incoherent.  We  are 
also  able  to  reproduce  the  RONAPP  detector  -esults  from  NORKSS  data, 
and  we  have  been  testing  a  NORESS  recipe  wi  h  72  coherent  and  12  Inco¬ 
herent  beams.  The  results  for  NORESS  with  t  ie  RONAPP  recipe  are 
exactly  the  same  as  for  RONAPP  itself,  when  using  same  detector  para¬ 
meters.  For  NORSAR  data,  the  results  are  better  than  the  old  NORSAR 
system,  as  coherent  beams  normally  have  20%  higher  signal-to-noise 
ratio  without  increasing  false  alarm  rate,  lost  false  detections  due 
to  spikes  in  data  are  eliminated  or  identif  ed  as  such.  Incoherent 
detection  results  for  NORSAR  are  equal  to  o:  better  than  in  the  old 
system.  These  better  results  are  obtained  even  though  new  delays 
within  subarrays  have  not  been  implemented.  In  Figure  4  some  of  the 
detection  results  are  listed. 
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So  far,  the  DP  system  has  be  :n  implemented  for  research  and  testing 
purposes.  The  code  contains  a  number  of  interactive  options,  and 
nothing  has  yet  been  done  t  >  reduce  overhead  for  data  administration. 
However,  the  beamforming  and  filtering  routines  are  coded  to  be  as 
fast  as  possible  with  FORTRAN.  In  addition  the  optimizer  option  is 
used  for  compilation.  Analyzing  the  time  consumed  by  the  different 
routines,  we  learn  that  beam  crming  is  the  most  time-consuming  part 
for  the  NORSAR  system.  Filtering  is  most  time-consuming  for  the  NORESS 
system.  To  run  the  NORSAR  recipe  for  180  coherent  beams  and  64  inco¬ 
herent  beams,  the  system  use  105X  of  an  IBM  4341/L01  or  40  %  of  an 
IBM  438L/P02.  The  relative  performance  of  IBM  4341/L01,  4341/M12 
and  4381/P02  for  adding  vectors  only  is  1.0,  1.5  and  2.6  respectively. 
The  requirement  for  doing  the  same  job  as  RONAPP  for  detection  only, 
is  about  1/3  of  the  NORSAR  processing.  An  early  version  of  the  program 
has  been  run  on  SUN/UNIX  system.  The  current  version  of  the  system  has 
been  especially  designed  to  avoid  possible  differences  between 
IBM/CMS  and  SUN/UNIX  systems. 

One  option  in  the  DP  system  that  we  have  found  very  useful,  is 
to  stop  the  system  at  some  s<  lected  time  point,  and  look  at  the 
current  data  buffers.  An  example  is  given  in  Fig.  IV. 5.  Any  of  the 
data  'channels'  may  be  displayed,  and  the  effect  of  varying  LTA  update 
parameters  may  be  visualized.  Likewise  different  STX  integration 
lengths  may  be  studied  to  set  the  capability  of  STX  as  a  measure  to 
pick  arrival  times. 

To  put  the  new  NORSAR  detection  system  into  operation  requires  only 
minor  changes  in  the  current  recording  system.  We  need  to  eliminate 
subarray  beamforming  and  filtering  on  the  M0DC0MP  computer,  and  we 
need  to  replace  detection  processing  on  the  IBM  side.  However,  the 
current  IBM  4341/L01  will  not  be  able  to  keep  up  with  real  time.  We 
may  solve  this  problem  In  twc  ways.  Either  by  trying  to  further  speed 
up  the  system  (some  overhead  may  be  reduced  by  making  a  copy  of  the 
program  which  is  able  to  process  NORSAR  data  only),  or  to  resample 
the  data  to  10  Hz  after  filtering  has  been  done.  The  latter  resampling 
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will  reduce  the  load  by  about  40  %.  We  will  look  further  into  these 
matters . 

The  DP  system  is  still  in  an  initial  version.  The  system  has  proved 
very  flexible  in  practice,  e.g.,  newly  started  research  on  using 
an  expanded  beam  set  for  the  NORESS  array  has  been  easily  adopted  by 
the  system.  However,  we  feel  that  there  is  still  room  for  improve¬ 
ments.  One  improvement  may  be  to  include  weighted  beamforming.  As  the 
program  may  easily  use  either  tape  or  NORSAR  online  disk,  we  may  com¬ 
pare  detector  results  for  any  period  of  time.  Inclusion  of  beamforming 
weights  or  other  detectors  can  be  easily  adopted  in  the  program. 

Better  performance  for  NORSAR  by  including  new  delays  on  the  subarray 
level  may  be  included  by  no  program  change,  as  this  is  input  through 
the  recipe.  Maybe  the  best  potential  for  improved  NORSAR  detector  per¬ 
formance  will  be  the  addition  of  high  frequency  digital  filters.  The 
computer  load  may,  however,  give  a  limitation  on  this,  at  least  in  the 
near  future. 


J .  F  /en 
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DATA  BUFFER 

LENGTH 

'.0.0 

SECONDS 

OLD  BUFFER 

LENGTH 

9 . 0 

SECONDS 

INPUT  THE  FOLLOWING 

STATIONS: 

01A01 

01A02 

01A03 

01B01 

01B02 

01  BO  3 

02B01 

02B02 

02  BO  3 

02C01 

02C02 

02C03 

03C01 

03C02 

03C03 

04C01 

04C02 

04C03 

06C01 

END 

* 

06C02 

06C03 

0 1A04SPZ 

01A05 

OlAOO 

01B04 

01B05 

01B00 

02B04 

02B05 

02B00 

02C04 

02C05 

02C00 

03C04 

03C05 

03C00 

04C04 

04C05 

04  COO 

06C04 

06C05SP2 

06C00 

SET  DETECTION  STATE  TRUE 


SET  REDUCTION  3.000 
SET  INITIAL  LTA  3.000 
SET  THRESHOLD  3.200 
SET  STX  LENGTH  0.500 
SET  STA  LENGTH  1.500 
SET  LTA  SIGMA  7.000 
SET  LTA  BETA  6.000 
SET  LTA  DELAY  9.000 


Seconds 

Seconds 

Seconds 

Seconds 


FILTER  FILNAM  FILTER  FLTTYP  BANDLW  BANDHG  IORD 


1 

FBI 

BU 

BP 

1.2 

3. 

.2 

3 

2 

FB2 

BU 

BP 

1.6 

3  . 

2 

3 

3 

FB3 

BU 

BP 

4 . 0 

6. 

,  0 

3 

FILTER 

ALL  WITH  FBI 

OUTPUT  xxxxxFBl 

FILTER 

ALL  WITH  FB2 

OUTPUT  XXXXXFB2 

DEFINE 

it 

NAO-AB 

AS  ALL 

xxxxxFBl 

USE 

NAO-AB 

USE 

i, 

BEAM-TABLE  NAO- 

451  TABLE 

* 

DETECT 

* 

ON  ALL 

DEFINE 

01A  AS 

01A01FB2 

01A02FB2 

01A03FB2 

01A04FB2 

01A05FB2 

0 1A00FB2 

END 

DEFINE 

01B  AS 

01B01FB2 

01B02FB2 

01B03FB2 

01B04FB2 

01B05FB2 

0 1 B00FB2 

END 

DEFINE 

02B  AS 

02B01FB2 

02B02FB2 

02B03FB2 

02B04FB2 

02B05FB2 

02B00FB2 

END 

DEFINE 

02C  AS 

02C01FB2 

02C02FB2 

02C03FB2 

02C04FB2 

02C05FB2 

02C00FB2 

END 

DEFINE 

03C  AS 

03C01FB2 

03C02FB2 

03C03FB2 

03C04FB2 

03C05FB2 

03C00FB? 

END 

DEFINE 

04C  AS 

04C01FB2 

04C02FB2 

04C03FB2 

04C04FB2 

04C05FB2 

04C00FB2 

END 

DEFINE 

06C  AS 

06C01FB2 

06C02FB2 

06C03FB2 

06C04FB2 

06C05FB2 

06C00FB2 

END 

*  NORSAR  Subarray  beam  #  1  for 

■  incoherent 

beams . 

MAKE  BEAM  01A01  OUTPUT  01A00SAB 
MAKE  BEAM  01B01  OUTPUT  01B00SAB 
MAKE  BEAM  02B01  OUTPUT  02B00SAB 
MAKE  BEAM  02C01  OUTPUT  O2C00SAB 
MAKE  BEAM  03C01  OUTPUT  03C00SAB 
MAKE  BEAM  04C01  OUTPUT  04C00SAB 
MAKE  BEAM  06C01  OUTPUT  06C00SAB 


USING 

01A 

VEL 

18.854 

AZI 

44 . 236 

USING 

01B 

VEL 

18.854 

AZI 

44.236 

USING 

02B 

VEL 

18.854 

AZI 

44 .236 

USING 

02C 

VEL 

18.854 

AZI 

44 .236 

USING 

03C 

VEL 

18.854 

AZI 

44.236 

USING 

04C 

VEL 

18.954 

AZI 

44.236 

USING 

06C 

VEL 

18.854 

AZI 

44.236 

MAKE  STX 
MAKE  STX 
MAKE  STX 
MAKE  STX 
MAKE  STX 


OF  01A00SAB  OUTPUT 
OF  OIBOOSAB  OUTPUT 
OF  02B00SAB  OUTPUT 
OF  02C00SAB  OUTPUT 
OF  03C00SAB  OUTPUT 


OlAOOSTXl 
01B00STX1 
02  B00STX1 
02CO0STX1 
03C00STX1 


Fig.  IV.  1  Extracts  from  the  NORSAR  detection  processor  "recipe". 

All  commands  except  the  complete  incoherent  recipe  are 
listed.  An  asterix  is  used  for  comments. 
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MAKE  STX  OF  04C00SAB  OUTPUT  04COOSTX1 
MAKE  STX  OF  06C00SAB  OUTPUT  06C00STX1 


DEFINE  NAO-INC1  AS  01A00STX1  OlBOOSTXl  02B00STX1  02C00STX1  03C00STX1 
04C00STX1  06C00STX1  END 


SET  DETECTION  STATE  TRUE 
SET  INITIAL  LTA  6.000 

SET  THRESHOLD  2 . 000 

SET  STX  LENGTH  0.500  Seconds 

SET  STA  LENGTH  1.500  Seconds 

SET  LTA  SIGMA  7.000 

SET  LTA  BETA  6.000 

SET  LTA  DELAY  9.000  Seconds 

* 

USE  NAO-IN  1 

* 

USE  BEAM  ■'•TABLE  NAO-271  TABLE  * 

* 

DETECT  INCOHERENT  ON  BEAM  1001 
DETECT  INCOHERENT  ON  BEAM  1002 
DETECT  INCOHERENT  ON  BEAM  1005 
DETECT  INCOHERENT  ON  BEAM  1006 
DETECT  INCOHERENT  ON  BEAM  1008 
DETECT  INCOHERENT  ON  BEAM  1010 
DETECT  INCOHERENT  ON  BEAM  1011 
DETECT  INCOHERENT  ON  BEAM  1012 
DETECT  INCOHERENT  ON  BEAM  1013 
DETECT  INCOHERENT  ON  BEAM  1015 
DETECT  INCOHERENT  ON  BEAM  1016 
DETECT  INCOHERENT  ON  BEAM  1018 
DETECT  INCOHERENT  ON  BEAM  1019 


Fig-  IV. 1  (page  2  of  2) 
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NAO-451  Bean  Table.  NORSAR  Coherent  beams. 


*  DELAYS  FOR  BEAM-TABLE  NAO-451  FOLLOW 

* 

STATIONS  01A01  01A02  01A03  01A04  01A05  01A00 

01B01  01 BO 2  01B03  01B04  01B05  01B00 

02B01  02B02  02B03  02B04  02B05  02BOO 

02C01  02C02  02C03  02C04  02C05  02C00 

O3C01  03C02  03C03  03C04  03C05  03C00 

04C01  04C02  04C03  04C04  04C05  O4C00 

06C01  06C02  06C03  06C04  06C0S  06C00  EN1) 

* 

SET  REFERENCE  02B00 


DELAYS  FOR  BEAM-TABLE  NAO-451 

* 


*  NAO  tf 

1 

VA 

10.890 

354 . 400 

GEO 

83 . 0 

-7.0 

C001 

42 

-1.665 

-1.268 

-1.132 

-1.382 

-1 .763 

-1.483 

-3.918 

-3.757 

-3.397 

-3.428 

-3.929 

-3 . 607 

-3.467 

-3.040 

-3.128 

-3.530 

-3.729 

-3.402 

-6. 124 

-5.674 

-5.603 

-5.963 

-6.415 

-5.977 

-5.658 

-5.223 

-5. 190 

-5.749 

-5.917 

-5.542 

-3 . 498 

-3.150 

-3.280 

-3.743 

-3.863 

-3 . 518 

-0.683 

-0.282 

0.000 

-0.254 

-0.765 

-0.441 

*  NAO  # 

2 

VA 

9.243 

35.800 

GEO 

73.0 

55.0 

C002 

42 

-0.817 

-0.480 

0.000 

-0.058 

-0.646 

-0.430 

-2 .458 

-2.592 

-2.171 

-1 . 868 

-2.228 

-2.138 

-4 . 004 

-3.536 

-3.  189 

-3.533 

-3 . 877 

-3.633 

-5.061 

-4 . 621 

-4.245 

-4.321 

-4 . 985 

-4.616 

-7.069 

-6.588 

-6. 154 

-6.541 

-6.977 

-6.670 

-6.305 

-5.812 

-5.587 

-5.935 

-6.321 

-6.009 

-2.159 

-1.779 

-1.272 

-1.177 

-1.719 

-1 . 602 

*  NAO  « 

3 

VA 

12 . 200 

324 .900 

GEO 

72.0  - 

74 . 0 

C003 

42 

-2.503 

-2.181 

-2.259 

-2 . 542 

-2.714 

-2.475 

-4 . 408 

-4 . 137 

-3.906 

-4 . 101 

-4.541 

-4.186 

-3.012 

-2.721 

-3.011 

-3.335 

-3 .405 

-3.123 

-5.769 

-5.435 

-5.541 

-5.951 

-6 . 1  2 

-5.816 

-4.207 

-3 .912 

-4.096 

-4 . 595 

-4.578 

-4 . 269 

-2. 140 

-1 . 961 

-2 . 240 

-2 .637 

-2 . 588 

-  2 .  318 

-0.394 

-0.090 

0 . 000 

-0. 367 

-0. 724 

-0 . 380 

Fig.  IV. 2  Extracts  trom  the  NORSAK  coherent  beam  table. 

The  first  3  beams  are  listed.  Beam  name  COo 1  -  CUUi 
corresponds  to  'old'  NORSAR  partition  1  beam  1  -  3. 
Delays  are  here  given  in  seconds  relative  to  the 
Instrument  with  latest  arrival.  The  program  will 
calculate  a  correction  time  for  each  beam  relative 
to  the  selected  reference  station. 
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DATA  BUFFER  LENGTH 
OLD  BUFFER  LENGTH 


60.0  SECONDS 
5.0  SECONDS 


INPUT 

THE 

FOLLOWING 

STATIONS: 

AOZ 

A1Z 

A2Z 

A3Z 

BIZ 

B2Z 

B3Z 

B4Z 

B5Z 

C1Z 

C2Z 

C3Z 

C4Z 

C5Z 

C6Z 

C7Z 

D1Z 

D2Z 

D3Z 

D4Z 

D5Z 

D6Z 

D7Z 

D8Z  D9Z  END 


USE  MAXAMPLITUDE  FACTOR  3.00  TO  MASK  DATA 

USE  NUMBER  OF  ITERATIONS  2  FOR  MAX  AMPLITUDE  MASKING 


SET  DETECTION  STATE  TRUE 
SET  REDUCTION 
SET  INITIAL  LTA 
SET  THRESHOLD 
SET  STX  LENGTH 
SET  STA  LENGTH 
SET  LTA  SIGMA 
SET  LTA  BETA 


3.000  Seconds 
20.000 
4 . 000 

0.250  Seconds 
1 . 000  Seconds 
6 . 000 
6.000 


SET  LTA 

DELAY 

5.000  Seconds 

FILTER 

FILNAM 

FILTER  FLTTYP  BANDLW 

BANDHG 

IORD 

1 

FB2 

BU 

BP 

1 .0 

3 . 0 

3 

2 

FB3 

BU 

BP 

1 .5 

3.5 

3 

3 

FB4 

BU 

BP 

2 . 0 

4.0 

3 

4 

FB5 

BU 

BP 

2.5 

4.5 

3 

5 

FB6 

BU 

BP 

3.0 

5.0 

3 

6 

BP4  1 

BU 

BP 

4 .0 

8.0 

3 

7 

BP4  2 

BU 

BP 

8.0 

16.0 

3 

8 

BP60 

BU 

BP 

1 .0 

2.0 

3 

9 

BP6J 

BU 

BP 

2.0 

3.0 

3 

FILTER 

AOZ 

WITH 

FB4 

OUTPUT 

A0FB4 

FILTER 

D1Z 

WITH 

FB4 

OUTPUT 

D1FB4 

FILTER 

D2Z 

WITH 

FB4 

OUTPUT 

D2FB4 

FILTER 

D3Z 

WITH 

FB4 

OUTPUT 

D3FB4 

FILTER 

D4Z 

WITH 

FB4 

OUTPUT 

D4FB4 

FILTER 

D5Z 

WITH 

FB4 

OUTPUT 

D5FB4 

FILTER 

D6Z 

WITH 

FB4 

OUTPUT 

D6FB4 

FILTER 

D7  Z 

WITH 

FB4 

OUTPUT 

D7FB4 

FILTER 

D8Z 

WITH 

FB4 

OUTPUT 

D8FB4 

FILTER 

* 

D9Z 

WITH 

FB4 

OUTPUT 

D9FB4 

* 

FILTER 

AOZ 

WITH 

BP60 

OUTPUT 

A0BP60 

FILTER 

C1Z 

WITH 

BP60 

OUTPUT 

C1BP60 

FILTER 

C2Z 

WITH 

BP60 

OUTPUT 

C2BP60 

FILTER 

C3Z 

WITH 

BP60 

OUTPUT 

C3BP60 

FILTER 

C4  2 

WITH 

BP60 

OUTPUT 

C4BP60 

FILTER 

C5Z 

WITH 

BP60 

OUTPUT 

C5BP60 

FILTER 

C6Z 

WITH 

BP60 

OUTPUT 

C6BP60 

FILTER 

C7Z 

WITH 

BP60 

OUTPUT 

C7BP60 

FILTER 

AOZ 

WITH 

BP6  1 

OUTPUT 

A0BP61 

FILTER 

C1Z 

WITH 

BP6 1 

OUTPUT 

C1BP61 

FILTER 

C2Z 

WITH 

BP6  1 

OUTPUT 

C2BP61 

FILTER 

C3Z 

WITH 

BP6 1 

OUTPUT 

C3BP6I 

FILTER 

C4Z 

WITH 

BP6 1 

OUTPUT 

C4BP61 

FILTER 

C5Z 

WITH 

BP6  1 

OUTPUT 

C5BP61 

FILTER 

C6Z 

WITH 

BP6 1 

OUTPUT 

C6PP61 

FILTER 

C7Z 

WITH 

BP61 

OUTPUT 

C7BP61 

Fig.  IV. 3  Complete  detection  "recipe"  for  emulating  RONAPP. 
(page  1  of  3) 
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DEFINE  TELEV  AS  AOZ 

C1Z  C2Z  C3Z  C4Z  C5Z  C6Z  C 7 Z 
!'1Z  D2Z  D3Z  D4Z  D5Z  D*Z  D7Z 
[•8Z  D9Z  END 


DEFINE  INTER  AS  AOZ 

HZ  B2Z  B3Z  B4Z  B5Z 
CIZ  C2Z  C3Z  C4Z  C5Z  C6Z  C7Z 
DIZ  D2Z  D3Z  D4Z  D5Z  D6Z  D7Z 
D8Z  D9Z  END 

* 

DEFINE  CRING  AS  AOZ 

BIZ  B2Z  B3Z  B4Z  B5Z 

CIZ  C2Z  C3Z  C4Z  C5Z  C6Z  C7Z  END 

* 

DEFINE  BRING  AS  AOZ  A1Z  A2Z  A3Z 

BIZ  B2Z  B3Z  B4Z  B5Z  END 


MAKE  BEAM  NRS1  OUTPUT  TELEV1  USING  TELEV  VEL  99999.9  AZI  0.0 

MAKE  BEAM  NRS2  OUTPUT  TELEV2  USING  TELEV  VEL  99999.9  AZI  0.0 

MAKE  BEAM  NRS3  OUTPUT  INTERS  USING  INTER  VEL  99999.9  AZI  0.0 

MAKE  BEAM  NRS4  OUTPUT  CRING4  USING  CRING  VEL  99999.9  AZI  0.0 

MAKE  BEAM  NRS5  OUTPUT  CRING5  USING  CRING  VEL  99999.9  AZI  0.0 

MAKE  BEAM  NRS6  OUTPUT  BRING6  USING  BRING  VEL  99999.9  AZI  0.0 

MAKE  BEAM  NRS7  OUTPUT  BRING7  USING  BRING  VEL  99999.9  AZI  0.0 

MAKE  BEAM  NRS8  OUTPUT  INTER8  USING  INTER  VEL  14.3  AZI  0.0 

MAKE  BEAM  NRS9  OUTPUT  INTER9  USING  INTER  VEL  14.3  AZI  90.0 

MAKE  BEAM  NRS10  OUTPUT  INTER10  USING  INTER  VEL  14.3  AZI  180.0 

MAKE  BEAM  NRS11  OUTPUT  INTER11  USING  INTER  VEL  14.3  AZI  15.0 

MAKE  BEAM  NRS12  OUTPUT  INTER 12  USING  INTER  VEL  14.3  AZI  75.0 
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V.  MAINTENANCE  ACTIVITIES 

y,  1 _ Activities  In  the  field  and  at  the  Maintenance  Center 

In  the  reporting  period  the  VMC  staff  have  been  Involved  in  many 
different  tasks  related  to  NORSAR,  NORESS,  and  the  second  array  in 
Finnmark.  The  activities  hav  comprised  corrective  and  preventive 
maintenance  for  NORSAR  and  NORESS,  including  the  hub  satellite  equip¬ 
ment.  Meetings  and  planning  related  to  the  second  array  (NORESS  II) 
and  continued  HF-experiments  in  connection  with  seismometer  tests  for 
the  HF-installation  were  als"  carried  out. 

At  the  NMC  laboratory,  major  tasks  have  been:  an  amplifier 
(3-component)  for  the  PDR-2  compuseis  Digital  Event  Recorder  built, 
repair  of  one  CCS-1  communications  system  equipment,  and  design  of  new 
anti-aliasing  filters  (B  Hz)  for  NORSAR. 

At  the  NDPC  NORSAR  communications  circuits  (comprising  modems  and 
lines)  have  been  regularly  monitored. 

The  NORSAR  subarray  electronics  have  been  checked  and  monitored  in 
accordance  with  established  .'.chedules. 

A  summary  of  the  activities  performed  during  the  period  is  given  in 
Table  V.l. 

For  completeness  we  list  activities  carried  out  for  both  systems, 
i.e.,  NORSAR  and  NORESS. 


Task 


Date 


Subarray/ 

area _ 

01A  Gain  SP  channel  01,  06  adjusted.  Mass  6  Oct 

position  all  LP  seism,  adjusted.  Test 
Generator  card  replaced. 

02B  Gain  SP  channel  01,  03,  04  and  Ob  adjusted.  7  Oct 

02B  (tele-  All  stations  visited  in  connection  with  in- 

metry)  spection  and  battery  water  refilling.  Cable 

channel  05  repaired. 

03C  Visit  in  connection  with  communication  ir-  l  Oct 

regularities,  had  been  unreliable  for  a 
while.  Several  attempts  were  made  to  find 
the  source,  also  NDPC  was  involved.  Data 
exchange  with  NDPC  OK  when  the  subarray 
was  left,  but  a  concrete  malfunction  was  not 
located. 

04C  Mass  position  all  LP  seism,  adjusted.  Ad-  15  Oct 

justed  free  period  (FP)  NS  seism.  +  gain. 

04C  Mass  position  all  LP  seism.,  NS  FP  and  gain  30  Oct 

adjusted. 

06C  Adjusted  SP  gain  ch  02,  04,  offset  ch  03  2  Oct 

and  mass  position  NS  seism. 

N0RESS  All  vertical  CS-13  seismometers  glued  to  Oct 

the  floor.  Water  leakage  in  the  cable  feed 
trough  -  tube  station  D5  stopped. 

Biasjd  Undenuten  station  visited  and  equipment  27,29  Oct 

repaired. 

NDPC  LP/SP-data  checked  dally.  Weekly  calibra-  Oct 

tion  of  SP/LP-instruments .  Offline  evalua¬ 
tion  of  03C  SP  channels  by  means  of  the 
CHANEV  SP-prograra.  Evaluation  job  ObC 
rejected  by  the  CHANEV  SP  program. 

02C  DC-offset  all  SP/ LP-channe Is  adjusted.  11,12  Nov 

Channel  gain  SP  02,  03,  0b  +  LP  vertical 
and  NS  adjusted.  Mass  position  and  FP  all 
LP  seismometers  plus  FP  vertical  and  NS 
seism,  adjusted. 
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Subarray/ 

area _ Task  _  Date 


SEtSNOR 

NDPC 

0  IB 

02B 

02C 

030 

04C 

06C 

NDPC 


NMC  participated  In  installation  of  seis-  18-20  Nov 
mometers  and  amplifiers 

Dally  check  of  SP/LP  data.  Weekly  calibra-  Nov 
tion  of  SP/LP  data.  SP-channels  all  subarrays 
evaluated  by  means  of  the  OFF-LINE  CHANEV  SP 
program. 


Adjustment  of  hP  channels  01,  04  and  03  3,12  Dec 

gain.  Adjustment  of  mass  position  (.ME)  on 
all  LP  seismometers. 

DC-offset  SP  channel  02,  03,  04  and  03  ad-  12  Dec 
justed.  Mass  position  all  LP  seismometers 
corrected.  Adjustment  of  free  period 
vertical  seismometers. 


Mass  position  IS  seismometer  adjusted.  1  Dec 

Cain  SP  channel  01,  DC-offset  all  LP  12  Dec 

channels  and  free  period  vertlcal/NS 
seismometers  adjusted. 

Modem  and  loop  control  logic  replaced.  11,12  Dec 

Adjusted  gain  SP  channel  03-06,  Vertical 
LP  channel  and  mass  position  (HP)  all  LP 
seismometers . 


Adjusted  mass  position  and  free  period  1  Dec 

all  LP  instruments. 

Daily  routine  checks  of  SP/LP  data.  Weekly  Dec 
calibration  of  3P/LP  seismometers.  A/D 
convertors  01A,  01B  evaluated  by  means  of 
collected  data  and  OFF-LINE  MtSN'i  program. 

SP  channel  01B  evaluated  after  data  collec¬ 
tion  and  the  OFF-LINE  CHANEV  SP  program. 
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Subarray/ 

area 

Task 

Date 

1987 

06C 

Preventive  maintenance  in  the  LP  vault 

13-15  Jan 
19,22  Jan 

NMC 

Built  a  3-component  amplifier  for  the 

PDR-2  recorder. 

Jan 

Much  time  spent  on  planning  in  connection 
with  the  second  array  in  Finnmark. 

NDPC 

Daily  check  of  SP/LP  data.  Weekly  calibra¬ 
tion  of  SP/LP  seismometers.  Besides,  OFF¬ 
LINE  program  M1SNO  run  on  subarrays  U2B, 

02C,  U3C  and  06C  in  connection  with  verify¬ 
ing  the  ability  of  the  A/D  conversion. 

Jan 

04C 

Modem  and  line  tests  carried  out  in  co¬ 
operation  with  NDPC. 

27  Feb 

NMC 

Two  PDR-2  compuseis  Digital  Event  Recorders 
and  1  CCS-1  compuseis  comm,  system  have  been 
repaired. 

Feb 

Planning  in  connection  with  the  second  array 
Finnmark  continued. 

in 

NDPC 

Daily  routine  checks.  Weekly  calibration  of 
SP/LF  instruments. 

Feb 

01A 

Communication  check  involving  modem  lines, 
etc.  Adjustment  of  gain  SP  channel  01,  U3 
and  06;  Mass  position  all  LP  seismometers, 
including  free  period  EW  seismometer. 

25,26  Mar 

NORESS 

Satellite  Tx  carrier  and  output  level 
adjusted. 

5  Mar 

NMC 

Daily  routine  check  of  SP/LP  data.  Weekly 
calibration  of  SP/LP  instruments. 

Mar 

Table  V.l  Activities  in  the  field  and  at  the  NORSAK  Maintenance 
Center,  Including  NDPC  activities  related  to  the 
NORSAR  array,  1  October  1986  -  31  March  1987. 
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y,2 _ Improvements  and  modifications  (NORSAR  arraj ) 

No  changes  since  last  reporting  period. 


V  .3 _ Array  status 

As  of  31  March  1987  the  following  channels  deviated  from  tolerances 


01A 

01 

8  Hz  filter 

02 

8  Hz  filter 

04 

30  dB  attenuation 

0  IB 

05 

04C 

01 

06C 

05 

Broadband  filter  installed. 

O.A.  Hansen 
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VI.  DOCUMENTATION  DEVELOPED 

Korhonen,  H. ,  S.  Pirhonen,  F.  Ringdal,  S.  Mykkeltveit,  T.  Kvaerna , 
P.W.  Larsen  and  R.  Paulsen:  The  Finesa  array  and  preliminary 
results  of  data  analysis,  University  of  Helsinki,  Inst,  of 
Seismology. 

Loughran,  L.B.  (ed.):  Semiannual  Technical  Summary,  1  April  -  30 
September  1986,  NORSAR  Sci.  Rep.  No.  1-86/87,  Kjeller,  Norway. 


L.B.  Loughran 
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VII.  SUMMARY  OF  TECHNICAL  RE PORTS/ PAPERS  PREPARED 

VT.t.l  Wide-band  slowness  estimation  using  a  small  aperture  seismic 
array 

In  order  to  evaluate  the  location  potential  of  the  NORESS  array,  the 
slowness  of  events  of  known  origins  have  been  estimated  applying  dif¬ 
ferent  approaches.  The  objective  of  this  investigation  is  to  develop 
methods  and  rules  that  can  be  applied  in  real-time  processing  of 
regional  array  data  and  to  give  an  idea  of  the  uncertainty  of  auto¬ 
matic  location  estimates. 

This  study  is  a  continuation  of  an  article  presented  by  Kvaerna  and 
Ringdal  (1986)  in  the  previous  semiannual  report.  In  that  work,  NORESS 
recordings  of  ten  quarry  blasts  from  a  dam  construction  site  in  south¬ 
western  Norway  (BlAsjo)  were  subjected  to  detailed  analysis  applying 
various  types  of  slowness  estimation  techniques.  Narrow-band  and 
wide-band  techniques  were  compared,  and  processing  in  fixed  frequency 
bands  were  tested  against  processing  in  the  band  with  maximum 
signal-to-noise  ratio  (SNR).  The  wide-band,  fixed  frequency  band 
approach  gave  by  far  the  most  stable  results.  The  maximum  spread 
around  the  mean  for  the  Pn-phase  was  only  1  degree.  By  that  approach 
we  were  also  able  to  separate  the  Sn  from  the  Lg  phases  on  the  basis 
of  estimated  apparent  velocity.  The  mean  azimuth  of  the  Lg-phases 
differed  only  1.4  degrees  from  the  theoretical  azimuth  and  the  maximum 
spread  around  the  mean  was  2  degrees. 

In  this  study  we  extend  the  data  base  by  events  from  two  other 
sites  with  known  location. 
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Semipalatinsk 

P-wave  recordings  from  sixteen  presumed  nuclear  explosions  at  Semipa- 
latinsk  have  been  analyzed  by  the  wide  band,  fixed  frequencies 
approach.  A  frequency  range  of  1.5-3. 5  Hz  was  selected.  This  is  for 
all  events  close  to  the  optimum  filter  band  for  SNR  gain  on  the 
array  beam. 

Site  coordinates:  50  N  79  E 

True  azimuth  from  NORESS:  75  degrees 
Distance  from  NORESS:  4200  km 

All  events,  except  one,  had  a  very  high  SNR.  As  shown  in  Kig.  VH.1.1, 
the  stability  of  the  Semipalatinsk  slowness  estimates  are  as  good 
as  in  the  Bias  jo  case.  The  mean  azimuth  is  78.10  degrees  and  the 
maximum  spread  around  the  mean  is  ±  1  degree. 

Estonia 

As  a  third  data  base  for  testing  of  the  location  algorithms,  NORESS 
recordings  of  ten  mining  explosions  in  the  Estonian  region  of  USSR 
were  chosen  (see  Fig.  Vlt.1.2).  These  events  were  located  by  the 
network  in  Finland,  and  reported  by  the  University  of  Helsinki. 

Site  coordinates:  59.5  N  28.5  E 

True  azimuth  from  NORESS:  92  degrees 
Distance  from  NORESS:  950  km 

In  the  5. 0-7.0  Hz  band,  the  Pn  phases  of  these  events  had  a  single 
channel  SNR  ranging  from  2.7  to  12.6,  which  is  much  lower  than  for  the 
B13sjs>  and  Semipalatinsk  events.  The  University  of  Helsinki  reported 
duration  magnitudes  (Md)  in  the  range  2.5  to  2.9. 
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Results  from  wide  band  slowness  estimation  of  the  Pn  phases  are  shown 
in  Fig.  Vtl.1.1.  They  were  all  processed  in  the  frequency  band  5. 0-7.0 
Hz.  The  mean  azimuth  is  81.38  degrees  and  the  maximum  spread  around 
the  mean  is  nine  degrees.  This  is  significantly  larger  than  for  the 
two  other  sites.  In  addition  the  mean  azimuth  has  a  bias  of  about  ten 
degrees  from  the  true  azimuth  of  92  degrees. 

Automatic  estimation  of  slowness 

It  is  difficult  to  predetermine  fixed  frequency  bands  tor  slowness 
estimation,  because  of  the  large  variation  in  signal  and  noise 
spectra,  see  Fig.  Vll.l.3b.  A  bad  choice  of  frequency  may  lead  to  pro¬ 
cessing  in  the  band  with  lowest  SNR  rather  than  best  SNK.  A  real-time 
detection  processor  will  use  several  filter  bands  and  array  sub- 
geometries  to  get  the  best  SNR  for  different  categories  ot  signals. 

The  detection  report  will  therefore  contain  information  about  which 
frequency  band  and  array  subgeometry  that  give  the  best  SNK .  By  using 
this  information  we  may  stabilize  the  slowness  estimates  In  various 
ways.  If  the  SNR  is  above  a  given  limit,  we  may  if  possible,  expand 
the  array  geometry  and  thereby  make  the  estimates  more  stable.  We  can 
also  stabilize  the  estimates  by  Increasing  the  processing  bandwidth  to 
higher  frequencies,  provided  the  SNR  is  sufflcently  high. 

We  have  tested  this  idea  by  processing  the  ten  Estonian  events 
with  the  same  array  suugeometry  and  filter  band  as  for  the  beam 
with  maximum  SNR.  The  results  are  shown  in  Fig.  Vll.l.la.  The  Pn 
and  Sn  esitmates  show  a  large  scatter  in  both  azimuth  and  velocity. 

In  case  of  Pn,  the  scatter  is  larger  than  that  shown  in  Fig.  VI  1.1.1, 
where  a  fixed  frequency  band  was  used.  In  comparison,  the  Lg  estimates 
are  quite  consistent.  The  mean  of  the  Lg  azimuths  is  92.98  degrees  with 
a  maximum  spread  around  the  mean  of  6  degrees.  We  were  not  able  to 
separate  the  Sn  and  Lg  phases  on  the  basis  of  velocity  at  in  the  Biasjd 
case,  probably  due  to  lower  SNR  and  a  narrower  usable  bandwidth. 
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Conclusions  and  recommendations 

This  study  has  shown  that  employing  a  wider  frequency  band  clearly 
tends  to  increase  the  stability  of  the  slowness  estimates,  provided 
the  signal-to-noise  ratio  is  adequate  over  the  band  of  interest.  The 
stability  was  found,  particularly  for  Pn,  to  be  remarkably  good  tor 
the  western  Norway  quarry  blasts  when  using  a  fixed  frequency  band  for 
all  ten  events.  The  same  conclusion  applies  to  the  teleseismic  p-phases 
from  Semipalatinsk. 

If  the  processing  band  for  slowness  estimation  is  selected  from 

the  detection  processor  information,  we  should  apply  the  estimated  Lg 

azimuth  for  location  purposes,  see  Fig.  VII. 1.3a. 

From  Fig.  VII. 1.3b,  we  can  see  that  the  Lg  spectra  show  a  more  uniform 
pattern  than  the  Pn  spectra.  The  best  frequency  band  for  detection  of 
Lg  phases  will  for  this  suite  of  events  be  1 .0-3.1)  Hz,  whereas  the 
optimum  Pn  filter  bands  will  range  from  2. 5-4.5  Hz  to  b.O-b.O  Hz. 

The  robustness  and  stability  of  the  Lg  azimuth  estimates  using  the 
detector  Information  are  related  to  the  fact  that  the  detection  pro¬ 
cessor  for  all  ten  events  reports  l.U-3.0  Hz  as  the  filter  band  with 
optimum  SNK.  Thus  the  detector-based  approach  is  in  this  case  equivalent 
to  the  wide-band  fixed  frequency  approach,  which  we  have  earlier  found 
to  produce  the  best  results. 


T.  Kveerna 
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Fig.  VI  1.1.1  Wide  band  slowness  estimates  of  the  P  phi.se  for  a  suite 
of  events  from  the  following  sites: 

i)  Sixteen  presumed  nuclear  explosions  from  Semipala- 
tinsk,  processed  in  the  1. 5-3-5  Hz  cand. 

ii)  Ten  explosions  from  a  mining  area  in  Estonia,  USSR, 
processed  in  the  5. 0-7.0  Hz  band. 

iii)  Ten  quarry  blasts  at  the  dam  construction  site  Bldsjo 
in  southwestern  Norway,  processed  in  the  5. 0-9.0 

Hz  band. 
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in  Che  Pn  spectra  of  the  upper  figure  with  two  peaks  at  about  i  Hz  and  6  Hz.  The  Lg  spectra 
lower  tigure  show  a  more  uniform  pattern  with  a  peak  at  about  2  Hz. 
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VIt.2 _ NORESS  noise  spectral  studies.  Noise  level  characteristics 

The  NORESS  noise  recording  system  has  continued  to  be  in  operation 
throughout  this  reporting  period.  The  system  is  documented  in  NORSAR 
Sci.  rep.  No  1-86/87,  where  also  results  from  the  first  half-year  of 
recording  are  presented.  The  system  is  used  to  obtain  statistics  on 
both  noise  suppression  by  beamforming,  and  the  actual  variations 
in  the  noise  level. 

The  beam  suppression  results  for  half  a  year  of  data  were  documented 
in  the  referenced  report.  The  continued  observation  of  beam 
suppression  has  reconfirmed  the  conclusions  made  in  that  report.  The 
results  for  the  different  subarray  configurations  are  consistent  and 
stable  over  time. 

We  will  therefore  not  report  any  further  on  the  beam  suppression  data, 
bu~  rather  concentrate  on  the  noise  level,  and  in  particular  extract 
features  important  for  signal  detection. 

An  initial  assessment  of  diurnal  and  weekly  variations  in  the  NORESS 
SP  noise  spectra  was  given  in  the  mentioned  report.  Observations  from 
more  than  one  year  of  data,  generally  confirm  the  previous  results, 
but  the  expanded  data  base  allows  us  to  describe  further  the  effect  of 
cultural  noise.  To  illustrate  time-dependent  varations  we  have  chosen 
to  use  uncorrected  data,  with  variations  in  power  measured  in  dB.  for 
illustration  of  the  actual  noise  level,  corrected  power  spectra  will 
also  be  presented. 

Pig.  VTT.2.1  displays  the  average  NORESS  SP2  noise  power  for  a  few 
selected  frequencies,  plotted  versus  time  for  the  four-week  period 
Monday,  15  December  1986,  through  Sunday,  11  January  1987.  A  vertical 
dotted  grid  line  is  plotted  at  00  GMT  for  every  24  hour  period,  in 
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this  figure,  days  which  are  not  normai  working  days  have  been  identi¬ 
fied  by  hatched  rectangles. 

As  may  be  seen  later  from  the  statistics  on  seasonal  variations,  the 
noise  level  during  the  Christmas  holidays  are  among  the  lowest  we  have 
observed  throughout  one  year  of  data  (for  frequencies  above  2  hz). 

We  note  that  although  industrial  plants  are  closed,  there  is  a 
cyclic  rise  in  noise  power  of  2  -3  dB  during  the  day,  with  a 
minimum  level  at  01  GMT  (local  time  02). 

For  days  with  normal  industrial  activity  the  diurnal  variation  may 
range  from  5  to  15  dB  for  frequencies  above  2  Hz.  However,  the 
nightime  minima  during  the  week  do  not  reach  the  minima  we  observe 
during  weekends.  These  observations  are  consistent  throughout  one  year 
of  data. 

These  data  represent  average  power  from  the  23  NORHSS  Sl’Z  instruments 
and  local  effects  from  wind  or  other  phenomena  which  may  disturb  one 
or  a  few  instruments,  will  be  smoothed  out  in  the  average  spectra. 

To  illustrate  the  seasonal  variations  we  have  chosen  to  use  obser¬ 
vations  made  at  particular  times  during  each  day.  For  nighttime  sta¬ 
tistics  we  have  thus  chosen  local  time  1)2  hours,  and  for  daytime 
observations  we  have  selected  local  time  13  on  normal  working  days 
only.  The  data  time  period  covered  is  31  March  1986  through  04  May 
1987.  Frequencies  selected  are  the  same  as  used  in  NOKSAK  Sci.rep. 

No.  1-86/87.  For  each  frequency  we  may  compare  observations  made  at  02 
and  13  local  time  to  see  the  effect  on  noise  level  from  cultural  acti¬ 
vities.  Results  are  presented  in  Figs.  VH.2  and  Vtt.2.3. 

Fig.  Vtt.2.2  shows  the  noise  levels  sampled  at  02  hours  local  time. 

In  the  figure,  each  vertical  dotted  grid  line  denotes  Monday  00  GMT  ot 
each  week.  There  are  7  observations  within  each  vertical  gridline.  We 
see  the  cyclic  variation  of  noise  level  reaching  a  minimum  every 


weekend,  and  a  pronounced  minimum  during  Christmas  holiday  and 
industry  vacation  (07  July  through  27  July).  Thus,  the  data  again  show 
that  the  variations  in  the  noise  level  above  2  Hz  are  to  a  large  ex¬ 
tent  due  to  cultural  effects.  Below  2  Hz  we  see  seasonal  variations 
that  are  not  reflected  at  the  higher  frequencies. 

Easter  holidays  (16  -  20  April, 1987  ),  are  also  characterized  by  rela¬ 
tively  low  noise.  However,  an  increase  in  noise  level  start  off  at  the 
end  of  the  holidays.  This  is  because  that  time  period  coincides 
with  the  snow  melting.  This  leads  to  greatly  increased  water  flow  in 
the  rivers  nearby  plus  transport  of  large  ice  blocks,  and  the 
corresponding  effects  on  the  seismic  noise  level  are  most  clearly  seen 
at  frequencies  above  1.5  Hz. 

The  snow  melting  period  in  Norway  normally  consists  of  two  phases: 
first  the  melting  of  snow  in  the  low-lands,  and  later  when  the  snow 
from  the  mountains  is  melting.  The  xlrst  period  strongly  affects  small 
local  rivers,  and  implies  a  sharp  increase  especially  in  high  fre¬ 
quency  noise.  The  second  period  causes  the  large  more  distant  rivers 
to  have  increased  water  flow,  and  the  effect  on  NORESS  is  mostly  con¬ 
fined  to  a  band  around  3  Hz. 

in  1986  these  periods  covered  the  interval  late  April  to  early  July. 
The  1987  snow  melting  started  about  20  April.  The  effect  of  water  flow 
in  the  rivers  was  also  documented  in  former  report. 

The  daytime  observations  are  shown  in  Fig.  Vll.2.3.  Many  of  the 
features  are  similar  to  those  observed  for  the  nighttime  data. 

Thus,  the  snow  melting  period  is  easily  detected.  Also,  the  holidays 
show  pronounced  minima,  with  Christmas  being  the  quietest  period. 
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Otherwise,  the  noise  level  at  daytime  is  consistently  above  that 
observed  at  night,  even  during  the  common  vacation  period. 

(07  July  -  27  July  1986). 

As  displayed  in  Fig.  VII. 2. 4,  data  from  the  NORESS  high  frequency  ele¬ 
ment  (here  named  HFZ),  show  the  same  features  with  respect  to  diurnal 
and  weekly  variations  as  do  the  standard  short  period  Instruments. 

We  may  infer  from  the  data  that  above  20.0  Hz  ,  the  industrial  noise 
is  not  as  pronounced  as  for  frequencies  in  the  range  2  -  20  hz. 
However,  as  this  data  is  from  from  one  single  instrument  as 
opposed  to  the  average  power  from  25  Instruments  for  NORESS  SPZ, 
the  data  show  much  more  variability.  Moreover,  the  instrument  is 
located  at  the  central  site,  closer  to  a  road,  and  high  frequency 
'signals'  from  local  activity  are  difficult  to  eliminate.  This  is  more 
easily  done  among  the  25  SPZ  instruments,  where  data  from  each  instru¬ 
ment  is  compared  to  the  average,  and  excluded  from  consideration  if 
interpreted  as  an  outlier. 

In  the  lower  part  of  Fig.  VI  1.2. 4,  we  have  displayed  uncorrected  power 
spcectra  for  the  high  frequency  element  observed  at  02  focal  time 
during  che  quiet  Christmas  period.  The  corrected  spectra  are  shown  in 
Fig.  VII. 2. 5.  The  peaks  in  the  spectra  are  due  to  electronic  inter¬ 
ference  from  fans  and  from  the  50  Hz  power  supply.  Fig.  VI  1.2.5  repre¬ 
sents  a  lower  bound  for  NORESS  noise  level. 

As  the  diurnal  and  weekly  variation  is  well  documented,  and  seasonal 
variations  are  mainly  tied  to  low-frequency  noise,  we  will  document 
the  actual  noise  level  at  the  NORESS  site  by  looking  at  ail  obser¬ 
vations  for  selected  time  points  throughout  the  last  two  half-year 
periods.  Again  we  look  at  local  time  02,  and  Vocal  time  Id  at  work¬ 
days.  The  spectra  are  corrected  for  system  response,  and  reported  as 
power  density  (nm**2/Hz).  Figs.  VII. 2. 6  -  VI 1.2. 8  show  the  results, 
both  for  NORESS  average  SPZ  and  the  High  frequency  element.  NORESS 
average  SPZ  data  Is  generally  consistent  witu  the  high  frequency  ele- 
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Bent,  but  small  variations  may  be  detected  due  to  the  averaging  of 
the  SPZ  instruments.  The  data  has  been  divided  into  the  two  half-years 
of  recording.  During  the  first  half-year,  the  high  frequency  data  was 
sampled  from  the  S-3  borehole  seismometer  at  a  depth  of  60  meter  (site 
F0).  During  the  second  half-year, the  data  is  from  the  GS-13  seismo¬ 
meter,  which  is  located  in  a  shallow  avult  (site  AO).  Comparing  the 
MEANZ  (average  SPZ  data)  with  HFZ  for  the  two  periods,  we  find  that 
the  MEANZ  and  HFZ  observations  are  close  to  equal  during  the  period 
when  GS-13  was  used. 

Comparing  the  upper  and  lower  parts  of  Fig.  VII. 2. 8,  we  see  some 
systematic  differences.  Thus,  for  frequencies  above  6.0  Hz,  the  sur¬ 
face  seismometer  shows  5  -  7  dB  higher  noise  level  than  the  borehole 
instrument . 


J.  Fyen 
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Fig.  VII. 2.1  Average  NORESS  SPZ  spectral  levels  plotted  on  an  hourly 
basis  for  the  four-week  period  15  december  1986  through 
11  January  1987.  Based  on  the  complete  noise  spectra 
sampled  each  hour,  the  plot  shows  the  power  versus  time 
for  frequencies  within  narrow  bands  (±  0.1  Hz).  Center 
frequency  for  each  band  is  printed  on  the  plot.  A  vertical 
dotted  grid  line  is  plotted  at  00  GMT  each  day  (01  local 
time).  Hatched  squares  Indicate  days  which  are  not  nor¬ 
mal  working  days. 
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Fig.  VII. 2. 2  Average  NORESS  SPZ  spectral  levels  observed  at  02  local 
time  only,  for  the  period  day  090,  198b  through  day  130, 
1987.  There  is  one  point  per  day  per  frequency  value,  and 
vertical  dotted  grid  lines  indicate  Monday  00  0MT  each 
week.  Frequencies  used  are  identified  by  the  center  fre¬ 
quency,  and  observations  that  are  within  ±  0.1  Hz  of 
center  frequency  are  used.  For  better  legibility,  the 
data  are  shown  in  two  separate  diagrams  (top  and  bottom) 
with  5  frequency  bands  in  each. 
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Fig.  VU.2.3  Average  NORESS  SFZ  spectral  levels  as  in  Fig.  VU.2.2,  but 
observed  at  13  local  time  on  working  days  only. 
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fig.  VU.2.5  NORESS  corrected  noise  power  spectra  based  on  data  from 
High  frequency  Element.  The  data  are  from  the  period  day 
359,  1986  through  day  004,  1987,  observed  at  0<i  local 
time  only. 
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Fig.  Vtt.2.8  NORESS  corrected  noise  power  spectra  as  in  Fig.  VU.i.^, 
but  for  local  time  13  and  working  days  only.  The  o<*ta  are 
from  the  period  day  091,  198b  through  day  273,  \98b 
(upper  part),  and  day  274,  1986  through  day  u9U,  1987 
(lower  part). 
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VII. 3 _ Initial  results  from  data  analysis  using  the  FINESA 

experimental  small  aperture  array 

This  contribution  describes  initial  research  results  obtained  in  con¬ 
junction  with  studies  of  data  recorded  by  the  FINESA  experimental  array 
in  Finland,  and  compares  these  results  to  those  previously  obtained  at 
NORESS. 


The  establisliment  of  FINESA  (Finnish  Experimental  Seismic  Array)  in 
late  1985  represents  a  joint  undertaking  between  the  University  of 
Helskinl  (Institute  of  Seismology),  Finland,  and  NTNF/NOKSAK,  Norway. 
All  data  and  results  from  the  array  are  openly  available  to  the  scien¬ 
tific  community. 

A  comprehensive  assessment  of  the  capabilities  of  FINESA  based  on  the 
first  year  of  operation  has  been  published  by  Korhonen  et  al  (1987). 

The  present  contribution  is  essentially  a  synopsis  of  that  paper, 
emphasizing  features  that  are  particularly  relevant  to  the  further 
development  of  the  small-array  concept  first  introduced  with  NORESS. 

The  reader  is  referred  to  the  referenced  paper  for  details  on  the  work 
described  in  the  following.  Copies  are  available,  upon  request,  from 
either  the  University  of  Helsinki  or  NTNF/NORSAK. 

Ar  ray  deployment 

The  location  of  FINESA  and  the  initial  array  geometry  are  shown  in 
Fig.  VII. 3.1.  The  array  is  situated  near  the  town  of  Heinola  about  100 
km  northeast  of  Helsinki,  and  the  distance  to  NORESS  in  Norway  is 
approximately  800  km.  FINESA  comprises  initially  a  total  of  10  SFE 
seismometers.  The  geometry  consists  of  a  3-instrument  ring  (A1-A3)  and 
a  6-instrument  ring  (B1-B6).  In  addition,  one  element  (Cl)  lies  outside 
these  two  rings.  The  maximum  intersensor  separation  is  about  1.5  km. 
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Tlie  A1  site  also  comprises  two  horizontal  SP  seismometers,  thus  pro¬ 
viding  a  3-corapouent  recording  system  within  the  array.  The  coor¬ 
dinates  of  the  individual  array  sites  are  given  in  Table  VI 1. 3.1. 

The  geology  of  the  siting  area  comprises  precambrian  gneissic  rock, 
with  little  or  no  sedimentary  overburden.  It  was  in  fact  possible  to 
place  all  seismometers  in  the  initial  deployment  directly  on 
outcropping  rock.  Topographically,  the  area  is  hilly,  with  the  maximum 
height  difference  between  sites  being  about  38  meters. 

All  12  seismometers  in  FTNESA  are  of  type  Geotech  S— 13.  Data  are 
transmitted  in  analog  form  by  cables  to  the  central  recording  site  A1 . 
The  recorder  is  a  Kinemetrics  DDS-1105  system,  which  comprises  12  bits 
linear  A/D  converters,  a  radio  receiver  for  time,  a  drum  recorder 
(type  Kinemetrics  VR1 )  and  a  magnetic  tape  transport. 

The  FINESA  system  has  been  in  operation  since  November  1485.  Most  ot 
the  time,  the  recording  system  has  been  operated  in  a  trigger  mode, 
using  a  built-in  voting  detector.  This  implies  recording  typically  two 
minutes  of  data  following  each  detection.  The  memory  hui  ter  ot  the 
Kinemetrics  system  is  somewhat  limited,  but  allows  3-4  seconds  ot 
noise  data  preceding  each  signal  onset  to  be  recorded  tor  all  chan¬ 
nels  . 

The  sampling  rate  is  selectable,  with  a  normal  setting  of  40  Hz.  The 
anti-alias  filter  prior  to  A/D  conversion  has  a  high  cutof f-f requency 
of  14.5  Hz.  A  calibration  signal  is  inserted  each  time  a  new  tape  has 
been  mounted,  thus  the  first  file  of  each  tape  is  a  calibration  file. 

Noi.8<;  studle£ 

For  array  sites,  it  is  important  to  establish  both  the  background 
noise  level  as  a  function  of  frequency  and  the  noise  correlation  pro¬ 
perties,  since  the  gain  that  can  be  obtained  by  beamforming  is  largely 
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dependent  on  the  structure  of  the  noise.  In  the  following  we  summarize 
our  findings  regarding  ground  motion  spectra  and  noise  correlation  by 
distance  for  the  FINESA  array,  and  make  comparisons  with  results  for 
the  NORESS  array. 

Ground  motion  spectra 

In  order  to  obtain  an  estimate  of  the  ambient  noise  level  at  the  site 
of  the  FINESA  array,  ground  motion  spectra  were  calculated  for  four 
time  intervals.  Two  intervals  were  from  November  1985,  one  during 
working  hours  and  one  during  night  time.  Two  intervals  from  October 
1986  were  chosen  in  the  same  way.  For  a  direct  comparison,  ground 
motion  spectra  for  data  recorded  at  the  NORtSS  array  were  calculated 
for  the  same  four  intervals.  Ail  noise  spectra  were  estimated  by  the 
indirect  covariance  method  described  in  Myktceltveit  et  al  (1985).  This 
method  gives  the  average  properties  of  each  noise  interval,  which  in 
our  computations  were  taken  to  be  100  seconds  long.  Each  noise  inter¬ 
val  was  visually  inspected  to  ensure  the  absence  of  definite  signal 
arrivals . 

Fig.  V1T.3.2  shows  ground  motion  spectra  for  two  of  the  four  intervals, 
for  both  FINESA  and  NORESS.  These  spectra,  and  other  spectra  shown  in 
this  report,  have  been  corrected  using  the  system  response  curves  for 
the  arrays. 

Our  studies  so  far  indicate  that  for  frequencies  up  to  2  Hz,  the  noise 

levels  at  FINESA  and  NORESS  are  comparable.  At  1  Hz  the  noise  level 

2 

ranges  between  2  and  20  nra  /Hz  for  all  noise  samples.  Between  U.5  and 
2.0  Hz  the  falloff  of  the  spectra  is  generally  steep,  ringing  from  23 
to  26  dB/octave  for  both  sites.  For  both  sites,  there  is  a  distinct 
change  in  the  spectral  slope  at  around  2.U  Hz.  For  frequencies  above 
2  Hz,  however,  the  shapes  of  spectra  at  FINESA  differ  from  those  at 
NORESS,  in  particular  during  night  time.  Between  2.0  and  3.0  Hz  the 
spectral  slope  has  decreased  to  7-10  dB/octave  for  FINESA  and  11-13 
dB/octave  for  NORESS.  The  four  FINESA  noise  samples  all  show  the  same 
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general  pattern  and  reach  a  minimum  level  of  42  to  46  dB  below  the 

2 

1  nm  /Hz  level  when  approaching  the  Nyquist  frequency  of  2U  Hz. 

Comparison  of  our  results  with  those  of  Bungum  et  al  (1465)  and 
Tarvainen  (1985)  confirm  that  FINESA  noise  level  at  high  frequencies 
is  relatively  high,  compared  to  some  other  sites  in  Fennoscandia ,  but 
otherwise  the  results  are  in  good  agreement. 

Noise  correlation  by  distance 

Both  for  NORESS  and  the  provisional  installations  preceding  its 
deployment  in  1984,  investigations  into  the  correlation  properties  of 
the  noise  field  played  an  important  role  in  design  and  evaluation 
efforts.  Curves  showing  intersensor  correlation  of  the  noise  versus 
sensor  separation  can  be  directly  utilized  in  predicting  the  signal- 
to-noise  ratio  gain  that  can  be  obtained  by  beamforming. 

We  have  computed  noise  correlation  curves  for  FINESA  in  the  same  way 
as  has  previously  been  done  for  NORESS  (see  Mykkeltveit  et  al,  1983; 
and  Bungum  et  al,  1985).  .Altogether  six  noise  intervals,  each 
comprising  110  seconds  of  data,  were  subjected  to  correlation  analysis 
in  six  frequency  bands.  For  each  of  the  frequency  bands,  the  six 
average  correlation  curves  were  plotted  on  top  of  each  other,  and 
"envelopes"  were  formed  that  contained  within  them  all  six  curves.  The 
results  for  four  of  these  bands  are  shown  in  Fig.  VII.  1. 3,  which  tor 
comparison  also  includes  envelopes  resulting  from  correlation  analysis 
of  noise  recorded  at  NORESS. 

All  FINESA  noise  correlation  curves  show  pronounced  negative  minima 
for  the  noise  correlation.  This  feature  has  been  shown  to  be  com¬ 
patible  with  a  model  of  isotropic  propagating  noise  near  the  Rayleigh 
wave  velocity  (Mykkeltveit  et  al,  1983).  Tt  has  also  been  demonstrated 
that  a  signal-to-noise  ratio  gain  beyond  the  standard  /N  gain  (N  being 
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the  number  of  sensors  in  the  array)  can  be  obtained  by  taking  advan¬ 
tage  of  the  negative  correiation  properties  of  the  noise,  and  the 
"deeper"  the  negative  minima,  the  larger  the  array  gain. 

The  FINESA  correlation  curves  appear  to  attain  their  minimum  values  at 
slightly  smaller  Interstation  distances,  compared  to  the  NOKESS 
curves  and  also  generally  reach  lower  negative  minima.  More  data  ana¬ 
lysis,  however,  is  needed  to  firmly  establish  this  conclusion.  The 
possible  Implication  is  that  an  optimum  array  at  the  FINESA  site  would 
have  a  slightly  smaller  aperture  than  an  optimum  array  at  the  NoRESS 
site,  and  that  such  an  array  might  have  potential  for  a  beamforming 
gain  even  exceeding  that  of  NORESS. 

Signal_ana  ly  si  s 

In  order  to  demonstrate  basic  array  signal  processing  techniques  and 
also  contribute  towards  an  assessment  of  the  potential  of  the  FINESA 
miniarray,  we  have  subjected  signals  from  several  regional  and  tele- 
seismic  events  to  detailed  analysis.  For  comparison,  NORESS  data  for 
two  of  the  regional  events  were  also  analyzed,  using  the  same 
algorithms.  Analysis  results  from  two  regional  events  recorded  at 
FINESA  are  shown  in  the  following;  otherwise,  reference  is  made  to 
Korhonen  et  al  (1987). 

Detailed  analysis  of  a  mining  explosion 

Fig.  VII. 3. 4  (top)  shows  the  FINESA  data  for  the  three-component  station 
at  site  Al  for  an  event  at  59.5’N,  28.5'E  on  December  27,  1983.  There 
are  two  distinct  P  phases;  the  first  arriving  P  phase  is  followed 
after  about  1.2  s  by  a  much  stronger  P  phase.  The  dominant  phase  in 
the  seismograms  is  the  Eg  (Sg)  phase.  We  also  see  clear  indications  of 
the  1  Hz  Rg  phase  in  the  Eg  coda.  In  the  bottom  part  ol  Fig.  VII.  1.4  we 
show  the  same  three  seismograms  low  pass  filtered  at  1  Hz,  and  the  Rg 
waves  stand  out  clearly. 
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Spectra  for  P,  Lg  and  Rg  for  this  event  are  shown  in  Fig.  VIT.3.5  along 
with  a  typical  noise  spectrum.  The  signal  spectra  are  computed  as 
described  in  Mykkeltveit  et  al  (1985)  and  are  based  on  data  window 
lengths  of  5  s  for  P  and  Lg  and  20  s  for  Kg.  We  see  that  all  phases 
have  considerable  energy  above  the  noise  level  throughout  the  fre¬ 
quency  band  1-20  Hz,  but  we  must  keep  in  mind  the  fact  that  the  Lg 
signal  window  used  contains  high  frequency  P  coda  and  likewise,  the  Rg 
window  contains  Lg  coda  arrivals. 

Phase  velocities  and  arrival  azimuths  were  estimated  using  the  wide¬ 
band  frequency  wavenumber  (f-k)  analysis  technique  described  by  Kvaarna 
and  Doornbos  (1986).  Fig.  VII. 3. 6  shows  frequency-wavenumber  spectra  for 
two  l  s  long  time  windows,  with  each  window  capturing  the  onset  of  a  P 
arrival.  Our  analysis  gives  phase  velocities  that  are  in  good 
agreement  with  results  from  refraction  profiling  experiments  in 
Finland  (e.g.,  Luosto  et  al,  1984),  for  the  appropriate  epicentral 
distance.  The  Lg  and  Rg  phases  were  also  subjected  to  wide-band  f-k 
analysis  and  results  are  shown  in  Fig.  Vll.3.7.  Again,  the  phase 
velocities  are  indeed  as  should  be  expected  for  these  phases.  We  see 
that  the  estimated  arrival  azimuths  are  within  2  degrees  of  the  true 
value  (calculated  from  the  network  location)  for  all  four  phases. 

NOKESS  data  for  a  three  component  station  for  this  raining  explosion 
(epicentral  distance  now  954  km)  are  shown  in  Fig.  Vtl.3.8  (top).  There 
are  three  distinct  phases:  Pn,  bn  and  Lg.  Spectra  for  these  phases  and 
the  preceding  noise  are  shown  in  the  bottom  part  of  Fig.  Vll.3.8.  We  see 
that  in  particular  Pn,  but  also  Sn,  has  energy  well  above  the  noise 
for  frequencies  up  to  15-2U  Hz.  In  comparing  FINESA  and  NORESS  P-wave 
spectra  for  this  event,  we  sec  that  there  are  similarities,  and  both 
spectra  exhibit  a  pronounced  peak  at  about  6  Hz. 
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Analysis  of  F1NESA  data  for  a  Bothnian  Bay  earthquake 

FINESA  data  for  an  earthquake  that  occurred  in  the  Bothnian  Bay  on 
October  27,  1986  (magnitude  =  2.6,  according  to  the  Helsinki  bulle¬ 
tin)  are  shown  in  Fig.  VII. 3. 9.  The  epicentral  distance  from  FINESA  is 
426  km.  The  phases  Pn  (very  weak,  just  above  the  noise  amplitude),  Pg, 
Sn  and  Lg  can  be  seen,  and  again,  Lg  is  the  overall  dominant  phase. 
Spectra  for  Pn,  Sn  and  Lg,  as  well  as  the  preceding  noise  are  shown  in 
Fig.  VII. 3. 9  (bottom  part).  We  observe  that  for  Pn  and  Sn,  the  SNR  is 
essentially  constant  over  the  frequency  ranges  8-13  Hz  and  4-15  Hz, 
respectively.  The  SNR  may  be  constant  over  an  even  larger  frequency 
range,  but  in  analyzing  these  data  we  are  limited  by  the  40  Hz 
sampling  rate  and  the  analog  filter  cut-off  frequency  of  14.5  Hz. 

Wide-band  f-k  spectra  were  computed  for  the  phases  Pn,  Sn  and  Lg  and 
results  are  shown  in  Fig.  VII. 3. 10.  We  see  that  again  the  resulting 
phase  velocities  are  in  agreement  with  the  expectations,  and  there  is 
also  an  indication  that  Sn  has  a  slightly  higher  phase  velooity  than 
Lg.  The  derived  azimuths  are  all  a  little  on  the  high  side  of  the 
"true"  value  of  323  degrees. 

IJvent  detection  studies 

In  the  initial  mode  of  operation,  the  FINESA  data  acquisition  system 
has  been  based  on  event  recording  using  the  built-in  detector  of  the 
Kinemetrics  DDS-1105.  A  detection  is  declared  whenever  the  short-term 
to  long-term  signal  ratio  exceeds  a  predefined  threshold  on  at  least 
two  out  of  three  pre-selected  channels. 

This  mode  of  operation  has  provided  a  considerable  amount  of  event 
recordings  for  research  purposes,  but  clearly  does  not  provide  a  good 
measure  of  the  inherent  detection  capability  of  the  array.  Therefore, 
we  have  also  studied  in  detail  continuous  recordings  over  a  shorter 
time  period,  using  off-line  processing  with  more  advanced  detection 
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algorithms.  The  results  from  this  experiment  are  described  in  the 
following. 


For  the  two-week  period  21  Oct  -  3  Nov  1986,  the  FINESA  array  was 
operated  in  the  continuous  recording  mode.  Because  of  the  frequent 
tape  changing  necessary,  the  standard  40  Hz  array  data  was  recorded 
only  during  daytime,  i.e.,  from  about  0400  GMT  to  about  1800  GMT.  For 
our  purposes,  this  was  quite  adequate  since  the  vast  majority  of 
reported  regional  events  (quarry  blasts)  occur  during  daytime. 

Detection  analysis  of  the  recorded  "continuous"  data  was  done  using 
the  R0NAPP  program  (Mykkeltveit  and  Bungum,  1984)  adapted  to  the 
FINESA  configuration.  Thresholds  were  set  equal  to  those  in  use  at 
N0RESS  and  this  resulted  in  a  relatively  large  number  of  false  alarms. 
Nevertheless,  actual  detections  could  be  confidently  Identified  by 
checking  signal  parameters  such  as  SNR,  frequency,  estimated  phase 
velocity  and  azimuth,  and  this  was  done  by  the  analyst  for  the  time 
period  2i  -31  Oct. 

The  detections  were  then  associated  to  the  Helsinki  regional  bulletin. 
Considering  only  those  bulletin  events  occurring  when  the  FINESA 
system  was  in  actual  operation,  we  found  that  out  of  68  reference 
events,  54  had  at  least  one  detected  phase  at  FINESA,  i.e.,  84  per 
cent.  This  could  be  compared  to  the  typical  NURESS  performance:  For 
one  year  (1985)  we  found  that  N0KESS  reported  at  least  one  phase  for 
28  per  cent  of  the  events  in  the  Helsinki  regional  bulletin  (996  out 
of  3521).  Thus,  the  advantage  of  being  close  to  most  epicenters  out¬ 
weighs  by  far  the  lower  SNR  gain  in  this  case. 

Both  the  P  and  S  (or  Lg)  phases  were  detected  for  the  ma  ;ority  of 
events.  In  addition,  several  Rg  detections  were  made  for  events  out  to 
about  250  km  eplcentral  distance. 
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In  Fig.  VII. 3. 11  the  azimuths  estimated  for  the  P  and  S  phases  are 
plotted  against  true  azimuth  (derived  from  the  Helsinki  bulletin).  The 
median  (50%)  error  is  about  5  degrees  for  both  P  and  S,  which  is  com¬ 
parable  to  initial  results  reported  for  NORESS  (Mykkeltveit  and  Ring- 
dal,  1981).  In  that  paper,  the  procedure  later  adopted  for  KONAPP  was 
used,  although  time  picks  were  made  by  the  analyst  rather  than  automa¬ 
tically. 

Fig.  VII. 3. 12  shows  estimated  phase  velocities  as  a  function  of  true 
azimuth  for  P  and  S  phases,  respectively.  The  line  at  6  km/s  is  the 
separation  line  currently  in  use  at  NORESS,  and  it  seems  to  provide 
quite  good  phase  discrimination  capabilities  between  primary  and 
secondary  phases  also  at  FINESA.  However,  the  separation  is  not  100 
per  cent,  and  there  is  a  fair  amount  of  scatter  in  the  estimated  velo¬ 
cities,  even  for  events  located  close  together.  We  consider  that  a 
somewhat  larger  array  geometry  will  be  necessary  to  obtain  more  stable 
phase  velocity  estimates,  although  broad-band  estimation  provides  some 
improvements  in  this  regard. 

Conclusions 

Noise  characteristics:  The  seismic  noise  level  at  the  FINESA  site  is 
slightly  higher  than  that  of  NORESS  at  high  frequencies  (2-15  Hz)  and 
comparable  to  that  of  NORESS  in  the  0-2  Hz  band.  This  is  consistent 
with  previous  observations  for  Fennoscandia .  No  particular  noise  sour¬ 
ces  have  been  identified  in  the  vicinity  of  the  array.  The  spatial 
noise  correlation  between  sensors  is  very  similar  to  that  observed  at 
NORESS;  in  particular  the  strong  negative  correlation  at  selected 
distance/frequency  combinations  is  present  also  at  FINESA,  and  appear 
to  be  even  more  pronounced  than  at  NORESS. 

Signal  characteristics :  The  signal  frequency  contents  at  FINESA  are 
comparable  to  those  observed  at  NORESS,  given  similar  epicentral 
distance.  Signal  correlation  across  the  array  is  high,  a^  expected, 
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except  at  very  high  frequencies.  A  noteworthy  feature  at  FINESA  is  the 
presence  of  strong,  coherent  Rg  waves  (i  second  period)  for  explosions 
out  to  more  than  200  km  distance,  in  contrast,  NORESS  does  not  observe 
Rg  beyond  50-60  km. 

Detection  capability:  The  FtNESA  detection  capability  is  excellent, 
both  for  regional  and  teleseismic  phases.  Our  initial  studies  have 
indicated  that  with  a  sensitive  KUNAPP-type  detector  the  array  is 
capable  of  detecting  more  than  80  per  cent  of  all  regional  events 
reported  in  the  University  of  Helsinki  Bulletin,  in  addition  to  a 
large  number  of  events  not  included  there.  The  detection  capability  is 
naturally  somewhat  lower  when  the  built-in  voting  detector  of  the 
present  Kinemetrics  system  is  applied. 

Phase  characterization:  The  phase  identification  capability  of  the 
initial  FtNESA  deployment  has  been  evaluated  using  both  single  fre¬ 
quency  and  broadband  f-k  analysis.  The  estimated  phase  velocities  are 
generally  reliable  indicators  of  phase  type  (primary  or  secondary), 
although  the  separation  is  not  100  per  cent.  The  resolution  is  not 
sufficient  to  separate  Sn  from  Lg  on  phase  velocity  alone,  a  result 
also  found  at  NORESS.  Azimuth  estimates  are  sufficiently  accurate  for 
successful  phase  association;  median  error  was  found  to  be  about  5 
degrees  both  for  P  and  S,  using  automatic,  single  frequency  f-k 
estimation. 

One  of  the  main  achievements  of  the  FINESA  study  has  been  to 
demonstrate  that  many  of  the  processing  techniques  that  have  contri¬ 
buted  to  the  outstanding  performance  of  NORESS  can  be  successfully 
applied  also  at  other  sites  in  Fennoscandla.  This  opens  up  for  the 


possibility  of  significantly  improved  capabilities  when  applying  data 
from  several  regional  arrays  in  a  future  multi-array  processing 
scheme. 
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Station 

Code 

Latitude 

ON) 

Longitude 

(,'E) 

ALNS 

61.4444389 

26.0793387 

ALZ 

61.4444439 

26.0793270 

ALEW 

61.4444462 

26.0793442 

A2 

61.4412337 

26.0736561 

A3 

61.4447211 

26.0715517 

B1 

61.4422706 

26.0844928 

B2 

61.4390893 

26.0984178 

B3 

61.4394023 

26.0705952 

B4 

61.4435798 

26.0650265 

B5 

61.4466901 

26.0710893 

B6 

61.4469378 

26.0808177 

Cl 

61.4384113 

26.0595485 

Accuracy  Internal: 
"  global: 


0. 1  m 

approx.  10  m 


Table  VII. 3.1  F1NESA  seismometer  coordinates. 


Altitude 

(.m) 


137.92 
137.8b 
137.92 
161.60 
133 .63 
164.98 
158.69 
175.57 

157.84 

142.68 

146.85 

157.69 
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Fig.  VII. 3. 2  FINESA  and  NORESS  ground  motion  spectra  for  day  322, 

1400  GMT  (top)  and  day  323,  0300  GMT  (bottom)  of  1985. 
The  FINESA  noise  spectra  are  given  by  dotted  lines, 
whereas  the  curves  for  NORESS  are  dashed-dotted. 


so 


IRSAR 


(985-36  t/(2:(6:48.4(9  59. SN  28. SE  MD  2 . 4 

□1ST.  264  KM  AZIMUTH  (49  DEG. 

_ _ _ _ _ 02/04/81  (2  ;Q4  ;Q  (  NO 


502. 


FA  1Z 

SPZ 


25  ■'» 


FAIN 

SPNS 


560. 


FA  IE 

SPEW 


SECONDS 

(985-56  (/ (2:  (6-.-i8 .4  (9  59. 5N  20. 5E  MO  2.4 

□  1ST.  264  KM  AZIMUTH  (4<J  OEG  . 

_ _  ...  .T'/m/Hi  ii;:(i:Q9  NORSAR 


Fiji*  VII. 3.4  FINtSA  data  tor  a  mining  explosion  (top)  and  the  same  dat 
low  pass  tittered  at  l  Hz  (bottom).  Distance  and  azimuth 
given  on  the  plots  are  according  to  the  epicenter  deter¬ 
mination  by  the  network.  ut  the  University  ot  Helsinki. 
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Spectra  tor  the  phases  F,  Lg  and  Kg  and  noise,  tor  dat 
shown  in  Fig.  VII. 3. d. 
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Fig.  VI  1.3.6  Wide-hand  f -k  spectra  tor  the  first  arriving  1'  phase  (top) 
and  flte  following  strong  P  phase  (bottom)  of  ttie  event  in 
Fig.  VII .  1.4. 
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-  84 


VT1.4 _ Multiple  scattering  by  topographic  relief 

Wavenumber  analysis  of  NORESS  array  data  has  revealed  sometimes  signi¬ 
ficant  frequency  dependent  effects.  Motivated  by  these  phenomena  we 
have  developed  a  method  to  account  for  multiple  scattering  due  to 
topography  of  internal  boundaries  or  of  the  free  surface.  The  method 
was  originally  developed  for  the  in  comparison  simple  problem  of 
electromagnetic  scattering  by  a  perfectly  conducting  tough  surface, 
and  it  is  based  on  the  so-called  extinction  theorem  te.g.,  brown, 
1985).  The  result  is  an  approximation  since  it  is  given  by  a  pertur¬ 
bation  series  whose  terms  are  obtained  recursively,  but  the  method  can 
be  applied  to  any  surface  tonography,  the  only  restriction  being 
implied  by  the  condition  of  convergence  of  the  series  solution.  The 
Born  approximation  is  simply  the  first  term  of  the  series;  it  is  thus 
straightforward  to  make  a  comparison  with  previous  single  scattering 
results.  Here  we  present  results  of  a  comparative  study  of  a  rough 
solid-liquid  boundary  where  explicit  single  scattering  results  are 
available  (Doornbos,  1978).  Topographic  effects  of  the  free  surface 
and  of  solid-solid  interfaces  will  be  considered  in  later 
contributions . 

The  solid-liquid  interface  forms  a  relatively  complicated  problem  due 
to  the  modified  boundary  conditions.  Let  the  deviation  of  the  inter¬ 
face  from  a  plane  be  given  by  z  =  f(x,y),  hence  the  reference  plane  is 
taken  to  be  the  z  =  0  plane  in  a  cartesian  coordinate  system.  The 
simplest  form  of  interfacial  displacement-traction  vector  needed  to 
determine  the  scattered  field  is 

d^  ■  (ux+,  Uy+,  vTu,  v^t v  /|_v|2)T  (1) 


t 
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Here  the  solid  is  taken  above  the  interface,  i.e.,  ux+  and  Uy+  are  the 
horizontal  displacement  components  in  the  solid.  The  vector  _v  is  nor¬ 
mal  to  the  interface: 


df  df 

(vx,vy,vz)  =  ( - , - ,  1) 

dx  dy 

hence  v^iv  is  the  normal  traction,  and  _v^u  is  a  constant  times  the 
normal  displacement  component.  The  solution  for  d  can  be  obtained 
recursively  in  wavenumber  space.  Let  I)(kx,ky)  be  the  Fourier  transform 
of  d(x,y),  and  let 


_D(kxfky)  31  £  D<n)(kx,ky) 

n=0 


(2) 


then  is  determined  from£(n-m),  m  *  l,...,n.  The  details  of  this 

solution  for  an  incident  plane  wave  are  giv  an  elsewhere  (manuscript  in 
preparation).  The  scattering  coefficients  are  similarly  deter¬ 

mined  fromJ)(n-m),  m  =  0,1, ...,n,  and  the  total  scattered  field  is 


JKkx»ky)  *  £  jl^n^(kx,ky)  (3) 

n-0 

The  components  of  are  just  the  plane  wave  reflection  and 

transmission  coefficients  for  a  plane  Interface,  jj(0)  is  the  Born 
approximation  representing  single  scattering,  and  B^n),  n  >  2,  include 
multiple  scattering.  From  the  components  of  1}  we  can  a.o.  calculate 
the  energy  flux  for  any  wave  type  j.  Per  unit  of  solid  angle: 
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2 

01 

EjO,*)  ~  ( - kx2-ky-)*  | Bj (kx ,ky ) | 2  (4) 

Vj2 


where 


(kx,ky)  **  —  (sin9  cos<t>,  sin8  sin  4>) 

vj 


and  Vj  is  the  wave  velocity. 

The  following  examples  are  for  a  velocity-density  structure: 

p+  =  5.56  g/cm3,  p-  =  9. ‘<3  g/cm3,  a+  =  13.64  Ijn/s,  a~  =  8.06  kra/s, 
P+  =  7.2  km/s,  =  0, 

a  space-stationary  roughness  function  f(x,y)  characterized  by  a 
Gaussian  autocorrelation: 


2  J 

r(x ,  y)  =  h  exp( - - —  ) 

a 

and  a  lateral  scale  length  =  20  km  lor  1  Hz  waves:  u>  =  2n  c/s. 

Fig.  VII. 4.1  illustrates  the  complete  series  solution  tor  incident  and 
scattered  wave  directions  which  are  typically  observed,  and  for  dit- 
ferent  topographic  height  h.  in  the  same  figure  we  show  the  phase  ot 
the  scattered  wave  due  to  a  single  bump  with  the  same  autocorrelation 
as  the  space-stationary  roughness  function.  The  first  order  phase 
solution  is  significantly  biased.  The  bias  turns  out  to  be  signiticant 
for  topographic  height  above  a  few  hundred  meters.  On  the  other  hand, 
the  first  order  approximation  predicts  the  scattered  energy  remarkably 
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well  even  for  rather  high  topography  (h  =  3  km).  Two  other  intuitively 
obvious  features  are:  (1)  The  low  order  solutions  are  relatively 
stable  for  wave  types  and  in  wave  directions  where  scattering  is  rela¬ 
tively  strong,  and  (2)  the  higher  order  terms  serve  to  enhance  energy 
relatively  far  from  the  specular  direction  and  to  reduce  energy  near 
the  specular  direction.  This  is  consistent  with  the  fact  that  multiple 
scattering  tends  to  widen  the  wavenumber  spectrum  of  the  scattered 
field.  The  difference  between  the  first  and  higher  order  solutions 
becomes  particularly  pronounced  in  directions  where  the  scattered 
field  is  relatively  small,  for  example,  near  a  zero  of  the  reflection 
coefficient.  Fig.  VII. 4. 2  illustrates  this  phenomenon  in  the  slowness 
range  above  3  s/deg. 


D.J.  Doornbos 


REFERENCES 

Brown,  G.S.  (1985):  A  comparison  of  approximate  theories  for  scat¬ 
tering  from  rough  surfaces.  Wave  Motion,  8,  155-205. 

Doornbos,  D.J.  (1978):  On  seismic  wave  scattering  by  a  rough  core¬ 
mantle  boundary.  Geophys.  J.K.  astr.  Soc.,  53,  643-662. 


A1S5  968  NORSAR  DETECTION  P 
COUNCIL  FOR  SCI  ENT 
L  8  LOUGH RAN  31  HR 
.9  IS  I  FI  ED  F08606-86-C-0004 


_  PROCESSING  SySTENIUj  R09AL  NORM] 

•UNCIL  FOR  SCIENTIFIC  AND  INDUSTRIAL  RESEARCH  K, 
5  LOUGHRAN  31  HAV  87  NORSAR-SCIENTTFIC-2-86-  87 


r 


-  90  - 
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Fig.  VII. 4. 2  Energy  flux  at  1  Hz  through  solid-liquid  interface  with 
the  same  specifications  as  in  Fig.  VlJ.x.l,  and  h  = 

2  km.  Scattered  P  in  the  liquid  from  incident  P  in  the 
liquid.  The  incident  and  scattered  ray  parameters  are 
equal.  (1):  Specular  direction.  (2):  The  vertical 
planes  through  the  incident  and  scattered  rays  differ  by 
44  degrees.  Solid  lines:  first  order  approximation. 
Dashed  lines:  complete  solution. 
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VII. 5 _ NORESS  capability  for  detection  and  location  of  mining 

tremors  In  the  Lubln  area  In  Poland 

The  small-aperture  NORESS  array  in  Norway  has  been  in  operation  since 
October  1984.  The  geometry  of  the  array  is  shown  in  Fig.  VII. 5.1.  The 
array  is  composed  of  25  short  period  vertical  seismometers  distributed 
in  four  concentric  rings  A.,  B,  C  and  D  around  the  central  site.  Since 
January  1985  data  from  NORESS  have  been  processed  in  real  time  at  the 
NORSAR  data  center  at  Kjeller.  The  data  used  in  the  detection  pro¬ 
cessing  comprise  25  short  period  channels,  sampled  at  a  rate  of  4U 
samples  per  second.  The  detection  and  location  algorithm  RONAPP, 
described  by  Mykkeltveit  and  Bungunt  (1984)  consists  of  digital  narrow- 
band  filtering,  beamforming,  linear  STA^LTA  detector,  frequency- 
wavenumber  analysis  and  association  of  regional  phases. 

Preliminary  evluation  of  the  NORESS  detection  and  location  capability 
of  regional  and  local  events  was  made  by  Mykkeltveit  et  al  (1985), 
Ringdal  (1986a,  b)  and  Mykkeltveit  (1986).  Rlngdal  (1986b)  concluded 
that  the  NORESS  detection  capability  appears  to  be  at  a  level  of 
magnitude  Ml  ”  2.5  or  better  up  to  a  distance  of  at  least  1500  km. 
Mykkeltveit  (1986)  studied  the  NORESS  automatic  performance  for 
seismic  events  in  Western  Norway  and  the  North  Sea,  situated  at 
distances  from  250  to  700  km  from  NORESS  and  within  a  magnitude  range 
from  1.4  to  3.1.  He  found  that  the  RONAPP  performance  for  events  at  a 
distance  of  400  km  from  NORESS  is  highly  satisfactory  for  events  with 
magnitude  2  and  greater. 

This  report  is  yet  another  attempt  to  evaluate  the  NORESS  real  time 
performance  related  to  small  events  Induced  by  mining  in  the  Lubin 
copper  basin  in  Poland.  The  Lubin  basin  is  situated  at  a  distance  of 
about  1060  km  from  NORESS,  in  a  tectonically  disturbed  area,  behind 
the  Tornquist-Teisseyre  tectonic  zone. 
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Data 

In  the  Lubin  basin  several  thousand  mining  tremors  are  recorded  and 
located  each  year,  with  local  magnitude  up  to  3.3  and  occasionally 
even  up  to  4.  All  seismic  events  are  recorded  by  three-dimensional 
underground  seismic  networks  operated  by  the  mining  industry.  They  are 
located  with  high  accuracy,  unobtainable  in  earthquake  seismology.  The 
most  accurate  parameters  are  the  epicenter  coordinates  with  an  error 
usually  not  larger  than  about  50-100  m.  The  depth  corresponds  in 
general  to  the  depth  of  mining  excavations.  The  least  accurate  para¬ 
meter  is  the  origin  time  which  is  not  determined  by  the  mine  networks, 
since  this  parameter  is  not  needed  for  mining  operations.  The  origin 
time  is  determined  for  larger  events  only,  in  most  cases  from  the 
records  of  KSP  station,  situated  at  a  distance  of  about  60  km  from  the 
mining  area.  Its  accuracy  is  not  better  than  about  0.2-0. 4  sec,  occa¬ 
sionally  worse.  Seismic  energy  rather  than  magnitude  of  each  located 
tremor  is  routinely  calculated  from  the  duration  of  the  seismograms 
recorded  by  mining  networks.  The  energy  can  be  readily  converted  into 
magnitude.  Additionally,  the  magnitudes  of  larger  events  are  computed 
and  verified  from  the  records  of  KSP  station. 

Between  January  1985  and  March  1986  three  hundred  and  eighteen  mining 
tremors  with  local  magnitude  from  2.0  to  3.4  have  occurred  in  the 
Lubin  basin.  A  complete  list  of  these  events  is  stored  in  the  NOKSAR 
filing  system.  The  NORESS  array  was  in  operation  during  the  occurrence 
of  265  events  (83  %  of  the  total  number,  Table  Vll.5.1),  and  it  was 
assumed  that  these  tremors  have  been  recorded,  regardless  of  their  SNK 
values . 

The  records  of  (a)  all  tremors  with  magnitude  Mg  >  2.5,  (b)  all  events 
detected  on-line  with  smaller  magnitude  and  additionally  (c)  25  smaller 
and  undetected  events  were  collected  on  separate  tapes  and  processed 
using  the  Interactive  computer  program  DISPAT  (Harri;  and  Kv*rna, 

1985).  Altogether  the  records  of  112  seismic  events  (Table  Vll.5.1) 
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were  processed  by  beamforming,  bandpass  filtering  and  spectral  analy¬ 
sis  of  P-  and  S-wave  groups.  As  a  result  of  the  previous  experience 
with  NORESS  records  of  regional  events  (e.g.,  Mykkeltveit  et  al,  1985) 
and  preliminary  experiments  with  the  records  of  Lubln  tremors,  a  fixed 
array  configuration  was  selected,  composed  of  the  center  instrument, 
C-rlng  and  D-ring.  The  processed  records  and  spectra  were  plotted  and 
then  compared  with  the  results  of  the  RONAPP  processing  on-line. 

Detection 

Three  data  sources  related  to  seismic  events  detected  by  NORESS  are 
available  at  NORSAR:  the  event  detection  list  and  standard  plots  of 
seismic  records  and  frequency  -  wavenumber  analysis,  produced  automa¬ 
tically  in  almost  real  time,  and  the  NORESS  event  bulletin  (available 
from  April  26,  1985)  which  is  based  on  the  detection  list,  verified 
and  processed  by  an  analyst.  These  sources  were  checked  against  our 
data  base  of  mining  tremors  from  the  Lubin  area.  It  was  found  that  68 
events  were  detected  automatically  by  RONAPP  (26  X  of  the  recorded 
number,  Table  Vlt.5.1).  The  number  of  recorded  and  automatically 
detected  events  as  a  function  of  magnitude  is  shown  in  Fig.  Vlt.5.2. 

In  57  cases  P  waves  were  detected  and  specified  as  such,  in  5  cases 
secondary  phases  (specified  by  RONAPP  as  Lg)  were  detected  as  first 
arrivals  and  in  4  cases  the  first  arrivals  were  unspecified.  In  2 
cases  the  detection  included  two  events  separated  in  time  by  one 
minute  (events  L  102  and  L  103)  and  by  half  a  minute  (events  L  337  and 
L  338).  Out  of  5  cases,  in  which  the  first  arrivals  were  specified  as 
secondary  phases,  in  3  cases  they  are  P  waves  (events  L  151,  L  161, 

L  251  with  magnitude  2.9,  2.7  and  2.5,  respectively).  Similarly,  the 
four  unspecified  first  arrivals  were  found  to  be  P  waves  (events 
L  267,  L  300,  1.  314  with  magnitude  from  2.2  to  2.6).  These  7  events 
are  not  listed  in  the  NORESS  event  bulletin.  The  arrival  azimuth  of 
recognized  P  waves  ranges  from  127*  to  187*  against  the  true  value  of 
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about  163’.  The  azimuth  of  P  waves  not  identified  by  RONAPP  as  such 
ranges  from  about  80’  to  300*. 

The  accuracy  of  automatic  determinations  of  P  arrival  times  was  tested 
by  comparison  with  the  arrival  times,  in  terms  of  travel  time  resi¬ 
duals,  read  directly  from  seismic  records  after  beam  forming  and  band¬ 
pass  filtering.  First  of  all,  the  average  apparent  velocity  of  P  waves 
between  NORESS  and  the  Lubln  area  was  calculated  from  46  selected 
records  with  clear  onsets  of  first  arrivals.  The  velocity  of  7.9 
km/sec  was  found,  and  this  value  was  accepted  for  the  calculation  of 
travel  time  residuals  (observed  -  calculated  travel  times).  The 
results  are  shown  in  Fig.  Vtl.5.3,  where  a  few  additional  observations 
from  the  RONAPP  processing  in  the  off-line  mode  of  undetected  events 
are  also  included.  The  residuals  calculated  from  direct  readings  range 
from  about  -1  to  1  sec,  while  those  corresponding  to  automatic  deter¬ 
mination  of  arrival  times  range  from  about  -1  to  3  sec.  The  largest 
difference  in  epicentral  distances  between  NORESS  and  the  Lubin  sour¬ 
ces  is  about  15  km.  Thus  the  arrival  times  picked  up  automatically  are 
on  the  average  l  sec  later  than  the  true  arrivals  and  often  they  are 
late  by  several  seconds.  This  delay  depends  strongly  on  the  shape  and 
corner  frequency  of  recorded  phases.  A  number  of  records  and  spectra 
are  used  below  to  illustrate  the  on-line  detection  and  location 
problems.  The  seismic  events  used  as  examples  are  listed  in  Table 
Vtl.5.2.  An  example  of  P-wave  records  from  event  L  311  with  the  arri¬ 
val  time  determined  automatically  about  1  sec  later  than  the  onset 
read  from  the  records  is  shown  in  Fig.  VII. 5. 4.  The  upper  trace  is  the 
NORESS  steered  beam  and  the  consecutive  traces  show  the  beam  filtered 
by  various  bandpass  filters  listed  in  Table  VII. 5. 3.  The  power  density 
spectra  for  the  beam  signal  of  6.5  sec  duration  and  beam  noise  of 
30  sec  duration  preceding  P  arrivals  are  also  shown  in  Fig.  VTT.5.4. 
The  spectra  were  obtained  by  a  prewhltenlng  of  the  data,  windowing  and 
transforming  the  auto-correlation  function  and  then  compensating  for 
the  prewhitening  (Harris  and  Kvaerna,  1986).  The  onset  picked  up  by 
RONAPP  corresponds  to  the  wave  packet  with  maximum  amplitudes  and 
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corner  frequency  of  about  3.5  Hz,  arriving  about  1  sec  after  the  first 
arrival  which  can  be  seen  at  lower  frequencies.  Another  example  of  P- 
wave  records  from  event  L  119  with  the  onset  determined  automatically 
about  2  sec  later  than  that  read  from  the  records  is  shown  in  Fig. 
Vtl.5.5.  Here  again,  the  automatic  determination  of  first  arrivals 
corresponds  to  high-frequency  wave  packets  with  maximum  amplitudes, 
although  the  real  first  arrivals  are  clearly  recorded. 

From  Fig.  VII. 5. 2  it  follows  that  no  clear  magnitude  threshold,  above 
which  all  events  are  automatically  detected,  can  be  established. 
Possibly  magnitude  2.7  approaches  such  a  threshold.  There  are  9  events 
with  Ml  “  2.6  and  10  events  with  Ml  >  2.b  which  were  not  detected.  Out 
of  the  9  events  with  Ml  *  2.6,  five  tremors  could  not  be  identified  by 
visual  inspection  of  the  traces  and  thus  are  not  expected  to  be  auto¬ 
matically  detected.  This  only  confirms  that  the  magnitude  value  of 
small  seismic  events  alone  is  not  a  sufficient  criterion  of  detection. 
The  noise  level,  dependent  on  the  season  of  a  year,  day  ot  a  week  and 
time  of  a  day  (e.g.,  Fyen,  1986),  is  a  major  factor  controlling  the 
SNR  values.  Four  events,  however,  with  Ml  “  2.6  are  clearly  recorded 
and  should  have  been  detected  (events  L  146,  L  169,  L  281,  L  286). 
Similarly,  out  of  10  undetected  events  with  Ml  >  2.6,  four  tremors 
with  magnitude  2.7  and  2.8  were  not  visually  identified,  but  5  events 
are  well  recorded  (events  L  112,  L  139,  L  175,  L  306,  L  353)  and  not 
detected.  The  sixth  event  L  104  with  Ml  “  3.0  occurred  i.5  minutes 
after  the  occurrence  of  major  tremor  L  103  with  Ml  *  3.4,  and  thi6 
possibly  prevented  its  separate  detection.  This  somewhat  erratic 
behavior  of  RONAPP  was  tested  by  standard  automatic  processing  in  the 
off-line  mode  of  23  selected  events,  well  recorded  by  N0RESS.  The 
assumption  was  that  possibly  RONAPP  is  not  always  in  operation  even 
when  NORESS  provides  undisturbed  recording  (RONAPP  operation  time  was 
not  available).  The  test  was  inconclusive,  out  of  23  trt mors  6  were 
detected  in  the  off-line  mode:  event  L,  169  with  Ml  *  2.6,  event  L  306 
with  Ml  •  2.9  and  4  events  with  smaller  magnitudes  (L  57,  L  60,  L  298 
and  L  302  which  was  also  located).  The  number  of  tremors  detected  by 
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RONAPP  in  the  off-line  mode  is  shown  in  brackets  in  Table  VIl.5.1.  Two 
larger  events  L  175  with  Ml  *  2.9  and  L  353  with  Ml  =  3.0  were 
detected  neither  in  the  on-line  nor  in  the  off-line  mode.  An  example 
of  such  a  case  is  given  in  Fig.  VII. 5. 6  where  the  P-wave  records  and 
spectra  from  event  L  112  with  Ml  =  2.7  are  shown. 

Occasionally,  however,  when  the  noise  level  is  exceptionally  low,  the 
RONAPP  package  performs  extremely  well  and  seismic  events  with  magni¬ 
tude  as  low  as  2.0  are  detected  (Fig.  Vll.5.2).  Such  a  case  pre¬ 
sented  in  Fig.  VII. 5. 7  where  the  P-wave  records  and  spectra  from  event 
1,  129  with  Ml  =  2.2  are  shown.  The  tremor  was  recorded  on  Saturday  in 
April  when  the  man-made  noise  with  a  maximum  peak  at  6  Hz  (e.g.,  Fyen 
1986)  was  absent.  This  is,  however,  another  example  of  the  automatic 
detection  of  first  arrivals  being  late  by  1.3  sec  from  the  true  onset 
time. 

Location 

Out  of  68  events  detected  on-line  by  RONAPP,  11  tremors  were  automati¬ 
cally  located  (16  per  cent  of  the  detected  number,  Table  VII. 5.1).  The 
number  of  detected  and  located  events  as  a  function  of  magnitude  Is 
shown  in  Fig.  Vlt.5.8,  and  the  list  of  located  events  is  given  in 
Table  VII. 5. A,  where  the  distance  between  the  NOKLSS  location  and  the 
true  epicenter  is  also  included.  The  essential  part  of  the  location 
procedure  is  the  association  of  P-  and  S-wave  arrivals  from  the  same 
event,  leading  directly  to  the  determination  of  an  epicenter,  since 
the  arrival  azimuth  is  already  known  from  f-k  analysis  (Mykkeltveit 
and  Bungum,  198A).  The  azimuth  determination  for  the  11  located  tre¬ 
mors  is  rather  consistent;  the  azimuth  ranges  from  1A9'  to  180'  for  P 
waves  and  from  158"  to  169'  for  S  waves  against  the  true  value  of 
163" •  P-  and  S-wave  arrival  times  from  12  located  tremors  (one  addi¬ 
tional  event  L  302  located  by  RONAPP  in  the  off-line  mode)  were  read 
diretly  from  the  records  and  compared  with  the  automatic  readings.  S-P 
intervals  from  automatic  detections  against  those  read  from  the 
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records  are  shown  in  Fig.  VII. 5. 9.  The  differences  between  those  two 
values  are  on  the  average  about  2  sec  and  occasionally  as  large  as  4 

sec,  which  corresponds  to  the  difference  in  epicentral  distance  of 

about  40-45  km. 

The  automatic  location  seems  to  be  the  most  unreliable  part  of  the 
RONAPP  processing  package,  and  the  problem  does  not  lie  in  the 
accuracy  of  the  location.  In  most  cases  the  secondary  phases  are  not 
recognized  by  the  system,  even  when  they  are  clearly  recorded,  and  the 

phase  association  part  of  the  processing  package  is  not  activated.  A 

number  of  such  cases  were  found  on  the  records  processed  with  DISPAT 
and  one  example  is  shown  in  Fig.  VII. 5. 10  where  well-recorded  S  waves 
from  event  L  82  and  their  spectrum  are  presented.  Naturally,  the 
secondary  phases  are  characterized  by  lower  corner  frequencies  than 
those  of  P  waves  and  only  the  records  filtered  with  lower  bandpass 
frequencies  display  distinct  S-wave  groups.  The  regional  event  loca¬ 
tion  algorithm  of  RONAPP  is  based  on  the  assumption  that  the  largest 
secondary  detected  phase  is  Lg  propagating  with  a  group  velocity  of 

3.5  km/sec.  This  assumption  is  valid  for  propagation  paths  confined  to 
continental  shield  type  structures,  like  Fennoscandia ,  but  is  not 
observed  for  propagation  paths  crossing  significant  tectonic 
disturbances.  Even  the  Tornquist-Teisseyre  tectonic  zone  in  Poland  is 
an  effective  barrier  to  the  propagation  of  Lg  waves  (Kvarna  and 
Mykkeltveit,  1986),  and  the  dominant  secondary  phases  recorded  by 
NORESS  from  the  Lubin  area  are  Sn  waves  propagating  with  a  velocity  of 

4.5  km/sec.  This  is  a  well  known  fact  and  further  studies  on  regional 
characteristics  of  secondary  phases  will  lead  to  a  vast  improvement  in 
the  event  location  of  regional  array  processing. 

The  epicenters  of  11  tremors  located  by  the  current  version  of  the 
NORESS  processing  package  are  shown  in  Fig.  VIT.5.11  (filled  circles). 
They  are  situated  at  a  distance  of  about  400  km  from  the  Lubin  area 
(Table  VII. 5. 4),  except  the  epicenter  of  event  L  258  whose  location  is 
based  on  a  misldentif led  secondary  arrival  by  35  sec.  These  epicenters 
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were  relocated  by  replacing  the  velocities  corresponding  to  Lg-Pn  time 
Intervals  by  those  corresponding  to  sn-Pn  intervals,  taken  from  the 
Jeff reys-Bullen  travel  time  tables.  The  results  are  also  shown  in  Fig. 
VII. 5. 11  (open  circles)  and  the  corrected  distances  between  the  new 
and  true  epicenters  are  listed  in  Table  Vll.5.4.  The  distances  are 
still  short  by  about  50-60  km  of  the  true  locations  as  a  result  of  the 
real  regional  wave  velocities  being  slightly  higher  than  those  from 
the  J-B  tables,  but  the  result  is  rather  impressive  considering  that 
it  is  obtained  from  a  single  small-aperture  array. 

The  NORESS  event  bulletin  gives  an  additional  indication  of  possible 
location  of  seismic  events  for  which  only  P  waves  are  detected.  The 
f-k  analysis,  performed  automatically  for  all  detected  P  waves,  pro¬ 
vides  the  arrival  azimuth  and  apparent  velocity,  which  in  turn  is  con¬ 
verted  into  the  epicentral  distance.  It  is  well  known  that  this 
approach  is  of  highly  limited  value  and  could  provide  only  very  rough 
indication  of  event  location,  and  should  not  be  discussed  here.  It  was 
found,  however,  that  the  records  of  P-wave  groups  display  charac¬ 
teristic  patterns  from  events  repeatedly  located  at  the  same  sites. 
This  fact  could  be  of  interest  for  an  expert  system  approach  to 
regional  array  processing. 

From  the  NORESS  bulletin  30  Lubin  events  were  found  to  be  located  from 
P  waves  alone.  Their  location  is  shown  in  Fig.  Vli.5.11  by  crosses;  9 
tremors  are  located  in  Roumania,  4  events  in  Bulgaria  and  Yugoslavia, 

4  events  in  Italy,  3  events  in  Czechoslovakia  and  4  tremors  are 
located  elsewhere.  Disregarding  the  last  four  events  which  are  ran¬ 
domly  distributed,  the  specific  locations  of  the  Lubin  tremors  seem  to 
be  related  to  characteristic  patterns  of  P-wave  records,  which  reflect 
the  source  properties.  This  problem  will  be  considered  elsewhere.  Here 
only  selected  examples  of  P-wave  patterns  are  given.  An  example  of 
characteristic  P-wave  records  and  spectrum  from  event  L  160  located  in 
Roumania  are  shown  in  Fig.  VI 1.5. 12.  The  largest  amplitudes  are 
observed  at  higher  frequencies  between  4.0  and  H.O  Hz,  about  2  sec 
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after  the  onset  (marked  by  the  first  arrow)  of  P  waves,  and  they  were 
recognized  by  RONAPP  as  P  first  arrivals  (marked  by  the  second  arrow). 
The  corner  frequency  is  5  Hz  (Fig.  VI1.5.i2c).  Another  example  of  P- 
wave  records  and  spectrum  from  event  L  337  located  in  Italy  is  shown 
in  Fig.  Vtl.5.13.  Here  the  largest  amplitudes  are  observed  at  lower 
frequencies  between  1.5  and  3.5  Hz  at  the  very  beginning  of  P-wave 
group  and  the  onset  of  P  waves  determined  automatically  is  very  close 
to  the  real  one.  The  corner  frequency  is  about  2.5  Hz  (Fig. 

VII. 5. 13c).  An  intermediate  case  is  presented  in  Fig.  VIT.5.14  where 
P-wave  records  and  spectrum  from  event  L  330  located  in  Bulgaria  are 
shown.  Here  the  maximum  amplitudes  are  observed  at  two  frequency 
bands,  between  1.5  and  3.5  Hz  and  between  4.0  and  8.0  Hz.  Two  corner 
frequencies  of  about  3  and  5  Hz  are  also  observed.  The  arrival  time  is 
about  1.7  s  later  than  the  onset  read  directly  from  the  records, 
though  the  first  arrivals  at  lower  frequencies  are  rather  distinct 
(Fig.  14b). 

Conclusions 

1.  A  data  base  of  318  mining  tremors,  with  magnitude  from  2.0  to 
3.4,  from  the  Lubln  area  in  Poland,  was  used  to  evaluate  the 
NORESS  real  time  performance  in  detection  and  location  of  small 
regional  seismic  events.  The  tremors  have  occurred  between  January 
1985  and  March  1986  at  a  distance  of  about  1060  km  from  NORESS. 

2.  The  array  was  in  operation  during  the  occurrence  of  265  Lubin  tre¬ 
mors  and  26  X  of  them  were  detected  on-line.  There  is  no  distinct 
magnitude  threshold  above  which  all  events  are  detected,  but 
magnitude  Ml  =  2.7  approaches  such  a  limit.  The  accuracy  of 
detected  arrival  times  and  azimuths  depends  strongly  on  the  shape 
and  corner  frequency  of  recorded  phases.  The  detection  system 
works  somewhat  erratically:  several  well-recorded  P  waves  were  not 
detected.  NORESS  is  obviously  a  unique  seismic  system  capable  of 
recording  so  small  events  at  a  distance  of  one  thousand  kilome¬ 
ters,  but  the  detection  system  certainly  needs  refinement. 
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Possibly  a  detector  operating  in  the  frequency  rather  than  in  the 
time  domain  would  be  more  suitable  for  the  Identification  of 
emergent  phases  and  those  recorded  with  low  SNR  (Gledhill,  1985). 

3.  The  location  algorithm  is  the  most  unreliable  part  of  the  pro¬ 
cessing  system  and  only  16  X  of  detected  events  were  located 
on-line.  Although  the  accuracy  of  location  is  rather  good  (within 
50-60  km  from  the  sourc"  area  after  repicing  Lg  by  Sn  velocities), 
it  can  be  further  improved  by  using  the  broad-band  f-k  estimation 
of  azimuths  and  velocities  instead  of  the  narrow-band  f-k  analysis 
presently  employed  by  the  KONAPP  package  (Kvaarna  and  Ringdai, 
1986).  In  most  cases,  however,  the  secondary  phases  are  not 
recognized,  even  when  tney  are  clearly  recorded,  and  the  phase 
association  part  of  RON.iPP  is  not  activated.  The  application  of 
horizontal  components  for  the  secondary  phase  detection  seems  to 
be  promising  (Kvserna,  1986),  but  it  is  not  clear  how  this  approach 
would  work  with  smaller  events.  A  rough  estimate  of  event  location 
from  P  waves  alone  could  be  improved  by  the  recognition  of  charac¬ 
teristic  patterns  of  P-wave  records,  reflecting  the  source  proper¬ 
ties,  from  events  repeatedly  located  at  the  specific  sites.  Such 
information  could  be  implemented  into  an  expert  system  approach  to 
regional  array  processing. 

4.  Obviously,  a  set-up  of  .mother  regional  array  of  the  NORESS  type 
would  vastly  improve  the  overall  capability  for  seismic  detec¬ 
tion  and  location  on-line.  At  the  present  stage,  however,  a  more 
vigorous  human  interaction  with  the  NORESS  automatic  system  Is 
badly  needed. 
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Magnitude 

ml 

Number  of  Events 

Reported 

Recorded 

Detected 

on-line 

Located 

on-line 

Processed 
with  DISPAT 

2.0 

54 

44 

3 

4 

2.1 

33 

28 

1 

2 

2.2 

41 

32 

2 

2 

2.3 

51 

42 

6 

14 

2  .4 

38 

32 

5  (3) 

(1) 

13 

2.5 

29 

25 

8  (1) 

15 

2.6 

22 

20 

11  (1) 

2 

20 

2.7 

20 

17 

13 

1 

17 

2.8 

11 

10 

6 

3 

10 

2.9 

10 

7 

5  (1) 

1 

7 

3.0 

7 

7 

5 

4 

7 

3.1 

1 

0 

0 

0 

0 

3.4 

1 

1 

1 

0 

1 

Total 

318 

265 

68  (74) 

11  (12)  112 

Per  cent 

83 

26  (28) 

16  (16) 

of 

of 

of 

reported 

recorded 

detected 

Table  VII. 5.1  Number  of  mining  tremors  from  the  Lubln  area  recorded 

and  detected  by  NORESS  from  January  1985  to  March  1986. 
The  number  of  events  detected  by  RONAPP  In  the  off-line 
mode  is  shown  in  brackets. 
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Center  Frequency 

Hz 

Filter 

Hz 

Band 

2.0 

1.0 

3.0 

2.5 

1.5 

- 

3.5 

3.0 

2.0 

4.0 

3.5 

2.5 

- 

4.5 

4.0 

3.0 

- 

5.0 

4.5 

3.5 

- 

5.5 

6.0 

4.0 

- 

8.0 

8.0 

4.0 

* 

11.0 

9.0 

6.0 

- 

12.0 

12.0 

8.0 

lb.O 

Table  Vtt.5.3  Bandpass  filters  used  tor  numerical 
filtration  of  beamed  records. 
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1  KILOMETER 


Fig.  VT 1.5.1  Geometry  of  the  NORESS  array.  The  four  rings  are  marked 

by  letters  A,  B,  C  and  D  and  the  four  three-component 
stations  are  marked  by  small  circles. 


Fig.  Vtl.5.2  Number  of  recorded  and  automatically  detected  events  as 
a  function  of  local  magnitude  Ml. 


TRAVEL  TIME  RESIDUALS  FOR  OFF-LINE  READINGS  (sec) 


Fig.  Vll.5.3  Travel  time  residuals  of  P  waves  calculated  for  automa¬ 
tic  determination  of  arrival  times  against  those  read 
from  the  records.  The  mean  values  are  shown  by  dashed 
lines. 
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Fig.  Vll.5.4  (a)  Example  of  F-wave  records  from  event  1311  (Table 

VII. 5. 2)  with  the  arrival  time  determined  automatically 
about  1  sec  later  than  the  onset  read  from  the  records 
(marked  by  arrows).  The  upper  trace  is  the  NURESS 
steered  beam  formed  from  the  center,  C-ring  and  D-ring 
vertical  instruments.  The  consecutive  traces  show  the 
beam  filtered  by  various  bandpass  filters  with  center 
frequencies  marked  at  the  records;  filter  bands  are 
listed  in  Table  Vlt.5.3.  The  time  scale  is  relative: 
the  absolute  P  arrival  time  corresponding  to  the  first 
arrow  is  given  in  Table  VII. 5. 2. 
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Fig.  VII. 5. 5  (a)  Example  of  P-wave  records  from  even.'  LI  1 9  with  the 

onset  determined  automatically  about  2  sec  later  than 
that  read  from  the  records  (marked  by  arrows). 

(b)  Power  spectra  for  the  beam  signal  of  7  sec  dura¬ 
tion  and  beam  noise. 
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Fig.  VII. 5. 6  (a)  Example  of  P-wave  records  from  event  LI  12  not 

detected  by  RONAPP.  (b)  Power  spectra  for  the  beam 
signal  of  8  sec  duration  and  beam  noise. 
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Fig.  VII. 5. 7  (a)  Example  of  P-wave  records  from  the  small  event  Ll.!9 

detected  automatically  as  a  result  of  unusually  low 
noise  level,  (b)  Power  spectra  for  the  beam  signal  ut 
b  sec  duration  and  beam  noise. 
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VII. 5. 8  Number  of  automatically  detected  and  located  events 
as  a  function  of  magnitude  ML. 
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Fig.  VTl.5.10  (a)  Example  of  well  recorded  S  waves,  event  Lt>2,  not 

detected  by  RONAFP.  (b)  Power  spectra  for  the  beam 
signal  of  9.5  sec  duration  and  beam  noise. 
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Fig.  VII. 5. 11  Geographical  distribution  of  seismic  events  from  the 

Lubln  area  (marked  by  a  square)  located  automatically  by 
NORESS  (marked  by  a  star).  The  epicenters  located  from 
P  and  secondary  phases  are  indicated  by  filled  circles 
and  their  positions  recalculated  by  replacing  Lg  by  Sn 
velocities  are  denoted  by  open  circles.  The  indications 
of  location  from  P  waves  alone  are  marked  by  crosses. 
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Fig.  VIZ. 5.12  (a)  Example  of  characteristic  P-wave  recoras,  event 
L160,  located  from  P  waves  only  in  Komania,  with  the 
largest  amplitudes  at  higher  frequencies  between  4.0  and 
8.0  Hz  (center  frequency  6.0  Hz),  (b)  P-wave  records 
scaled  to  the  maximum  amplitude  at  the  center  frequency 
of  6.0  Hz. 
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Fig.  VII. 5. 13  (a)  Example  of  characteristic  P-wave  records,  event 

L337,  located  from  P  waves  in  Italy,  with  the  largest 
amplitudes  at  lower  frequencies  between  1.5  and  3.5  Hz 
(center  frequency  of  2.5  Hz),  (b)  P-wave  records  scaled  to 
the  maximum  amplitude  at  the  center  frequoncy  of  2.5  Hz. 
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Fig.  VU.5  .14  (a)  Example  of  characteristic  P-wave  records,  event 

L330,  located  from  P  waves  in  bulgaria,  with  the  largest 
amplitudes  at  two  frequencies,  between  1.3  and  3.3  Hz 
(center  frequency  2.3  Hz)  and  between  4.0  and  b.O  Hz 
(center  frequency  6.0  Hz),  (b)  P-wave  records  scaled  to 
the  maximum  amplitude  at  the  center  frequency  of  6.U  Hz. 
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